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conventional  survey  methods  is  elevated  to  absolute  accuracy 
when  redundant  observations  are  made  and  adjusted  using  the 
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When  redundant  observations  are  available,  hydrographic 
positioning  problems  require  the  application  of  a  data  adjust¬ 
ment  method  so  that  all  information  may  be  used  for  obtaining 
the  most  reliable  ^fix* i  One  of  the  oldest  and  best  engineering 
techniques  developed  for  the  purpose  is  based  on  the  least 
squares  principle.  The  theoretical  background  is  provided  to 
explain  that  principle  and  the  technique  for  its  application. 
Also,  the  analytical  solutions,  and  respective  computer  programs 
implementing  them,  are  developed  for  the  following  hydrographic 
positioning  methods: 

>  a)  fix  by  N  azimuths;- 
-b)  f^Lx  by  N  sextant  angles'  &  L 
- c)  fix  by  two  range  distances  and  one  azimuth. 

For  each  method,  an  illustrative  application  of  the  respective 
computer  program  is  presented. 

The  least  squares  adjustment  method  not  only  yields  the 
most  likely  values  for  the  fix  coordinates  but  also  statistically 
quantifies  position  accuracy.  Relative  accuracy  achieved  with 
conventional  survey  methods  is  elevated  to  absolute  accuracy 
when  redundant  observations  are  made  and  adjusted  using  the 
method  of  least  squares. 
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In  hydrographic  surveying,  the  determination  of  position 
is  as  important  as  the  measurement  of  depth.  Conventional 
survey  methods  rely  primarily  on  two  lines  of  position  (LOP) 
to  establish  a  fix.  These  LOP's  are  obtained  by  measuring 
angles  and  distances  directly.  Alternately,  electronic 
positioning  systems  are  used  to  establish  a  pattern  of  LOP's 
(an  electronic  lattice)  based  on  the  propagation  of  electro¬ 
magnetic  energy. 

Until  recently  it  has  been  logistically  unfeasable  to 
obtain  redundant  observations  in  hydrographic  surveying. 
However,  with  the  advent  of  computers  and  miniaturized 
electronic  positioning  systems,  redundant  observations  are 
now  being  made  in  order  to  increase  fix  accuracy  and  prevent 
delays  due  to  equipment  malfunction. 

Mathematical  adjustment  methods  must  be  applied  to  the 
redundant  data  sets  in  order  to  maximize  the  accuracy  of 
each  fix.  One  such  adjustment  method  is  based  on  the 
principle  of  least  squares.  It  assumes  that  blunders  and 
systematic  errors  have  been  removed  from  the  data  so  that 
only  random  errors  remain.  This  method  yields  not  only  the 
best  estimate  of  position  for  a  given  set  of  redundant  LOP's, 
but  also  assesses  the  absolute  accuracy  associated  with  each 
fix  determination. 
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A.  INTRODUCTION 


In  general,  the  redundant  observations  of  any  variables 
in  a  physical  system  (such  as  in  hydrographic  position 
determination)  do  not  precisely  satisfy  the  mathematical 
model  developed  to  represent  that  system.  However,  the 
derivation  of  every  mathematical  model  is  based  on  the 
assumption  that  the  true  values  of  the  variables  will  satisfy 
the  model.  The  difference  between  the  true  and  observed 
value  for  any  physical  variable  is  called  the  residual; 

residual  =  t rue  value-  observed  value,  (i-l) 

In  making  physical  measurements,  true  values  can  never 
be  determined.  Considering  the  observed  values  as  values 
assumed  by  random  variables  following  normal  distributions, 
every  true  value  can  be  represented  as  the  mean  of  a  random 
variable.  Therefore,  eq.  (1-1)  can  be  rewritten  as 

residua  l »  mean  of  random  varia  ble-  ohs.  Value  .  ( I-Z) 

The  least  squares  principle  establishes  a  criterion 
for  obtaining  the  best  estimates  of  the  true  values.  It 
states  that  the  true  values  will  be  such  that  the  sum  of 
squared  residuals  is  a  minimum.  For  a  further  discussion  of 
this  principle,  see  Appendix  A. 
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B.  LEAST  SQUARES  PRINCIPLE  FOR  UNWEIGHTED  OBSERVATIONS 
Measuring  different  parameters  of  a  mathematical  or 
functional  model,  we  associate  with  each  parameter  a  random 
variable, .  Designating  by  the  value  assumed  by  the 
random  variable  (the  observed  value)  and  by  pi  its  mean 
(the  adjusted  value) ,  the  residual  is  given  by 


For  Jl  observed  parameters,  the  least  squares  fundamental 
condition  is  expressed  as 

H  „ 

Z  /Z  i  ,1  /t 

Yl  s  V,  ¥  *1  +  *  •  .  +  Yj,  S  minimum 

i  •/ 


or,  in  matrix  form 


\J  \J  ~  minimum 

where  V  *  ^  V,  .. . 


(1-4) 


C.  LEAST  SQUARES  PRINCIPLE  FOR  WEIGHTED  OBSERVATIONS 

If  the  fl  observations  are  unequally  weighted,  then  the 


least  squares  fundamental  condition  is  expressed  as 

t* 

->  A 

'L 

i » i 


n 

Vi  -  i  \(  +•••  +  Yj. 
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or,  in  matrix  form 


V  WV  =  minimum  (  1-5) 

where  W is  the  ff*h  weight  matrix.  See  Appendix  B  for  a 
more  complete  discussion  on  the  concept  of  weighted 
observations . 

D.  OBSERVATION  EQUATIONS 

In  the  expression  for  the  residual,  ^  is  a  known  value 
(the  observed  value)  and  represents  (from  a  deterministic 
point  of  view)  the  true  value,  thus  satisfying  the  relation¬ 
ship  between  the  variables  as  expressed  in  the  functional 
model.  The  model  must  define  an  analytical  expression 
relating  the  unknown  values  with  the  known  ones.  In  general, 

Jl l  may  be  expressed  as  a  function  of  the  unknowns; 

s  j’i  (*i  •*  • 

where  X,  are  the  unknowns.  Therefore,  eq.  (1-3) 

can  be  rewritten  as 

Vi  *  ft  (*.  *x  •••  "m  )-  Yc.  (  I-  Q 

The  above  expression  is  referred  to  as  an  observation  equation. 


13 


If  is  a  linear  function,  the  observation  equation 
may  be  written  as 


Vl  *  3 Lo  +  lit  *i  +  •  •  *  +  X m  ~  jfi  L  ?) 


where  3io  ;  j  3^  are  coefficients.  The  least 

squares  method  does  not  require  that  the  observation 
equations  be  expressed  in  linear  form.  However,  the  compu- 
tations  to  determine  the  values  of  the  unknowns  are  greatly 
simplified  if  the  observation  equations  are  linearized. 

If  fa  is  nonlinear,  a  Taylor's  series  expansion  may  be 
applied  to  linearize  the  function.  Since  it  is  not 
practical  to  work  with  all  the  terms  of  the  expansion,  only 
the  zero  and  first  order  terras  are  used.  Thus,  the  linear¬ 
ized  form  is  an  approximate  analytical  expression  for  fa ; 


fl‘  k+  ••••***-"  )  ft L  ■ 

Since  the  function  fa  and  its  partial  derivatives  may  be 
evaluated  given  approximate  "initial  values"  for  the  unknowns, 
the  observation  equations  can  be  expressed  as  linear  functions 
of  the  increments 


+  3  tin 


■ft  «  aio  +  ^x,  ¥ 


4X,n 


where 


3ic  *  ji  ^ 


and 


3Ct  »  tli 


Therefore,  the  residual  Vt  will  be  stated  as 


3to  +  ax',  +. . 


+  3  C  m  A  X  m 


It  must  be  emphasized  that,  using  the  approximate 
expression  for  ,  the  least  squares  method  will  yield 
adjustments  (  A*,  •  •»  AVm  )  which  must  then  be  applied 

to  the  "initial  approximations" . 

Therefore,  an  iterative  solution  is  required  to  solve 
for  the  final  values  of  the  unknowns.  The  first  adjusted 
results  are  used  as  the  new  initial  values,  and  the  obser¬ 
vation  equations  must  be  formulated  again.  This  process  is 
continued  until  the  increments  become  vanishingly  small  or, 
from  a  practical  point  of  view,  converge  to  within  a 
specified  tolerance. 


E.  LEAST  SQUARES  ADJUSTMENT  METHOD 

Considering  eq.  (1-7) ,  and  combining  the  constant  terms 
a  new  expression  for  the  observation  equation  is  obtained 


L5 


Then,  the  n  observation  equations  can  be  presented  as  the 
following  system  of  ft  equations  with  tn  unknowns,  where 
n>m  for  the  case  of  redundant  observations: 


Y  “  X  i  +  3im  "*  ^  / 


V*|  «3f»i  X j  +  •  •  .  +■  Smk  -  in 


fl-  10 


or,  in  matrix  notation,  as 


V*  AX  -  L  (i-n 


where 


4. 


3»(  •  '* 


Equations  (1-10)  and  (1-11)  express  the  general  form  of  the 
observation  equations.  By  imposing  the  least  squares  prin¬ 
ciple  that  VTVV  V* minimum,  a  set  of  equations  are  obtained 
which  can  be  solved  to  find  the  best  estimate  of  the  unknown 
values.  These  expressions  are  known  as  the  normal  equations. 
For  the  observation  equations  as  expressed  in  eq.  (1-10) , 
they  form  a  set  of  m  equations  and  m  unknowns: 

[«U£  dil  ]#*♦•.  •  +  [**£  ]  Xm  li.  ]  *  O 

. _ .  (M2) 

|w>;3in*  •  •  +  [u K  aim**  ]  Xm-  ^  cue  aim  U  j  s  o 

where  the  brackets  have  the  usual  meaning  of  sum  (  L  * 

In  matrix  notation,  the  normal  equations  are  expressed  as 

(ArWA)  X  -  i^TW  L.  (1-13) 

The  normal  equations  are  used  to  solve  for  the  values  of  X J 

X  *  (aVa)'  (4tW  L)  (1-14) 


where  X  is  the  vector  whose  elements  are  the  adjusted  values 
for  the  unknowns.  For  a  more  complete  development  of  the 
normal  equations  see  Appendix  C  and  Appendix  D. 


F.  PRECISION  OF  OBSERVATIONS 

When  observing  an  unknown  variable  a  finite  number 
of  times  t  n  ,  the  value  of  can  be  estimated  by  comput¬ 
ing  a  sample  standard  deviation ,  S  ,  according  to  the 
following  formula: 


where  Y^  *,h)  represents  the  observed  values  and 
X  the  average  value,  for  a  set  of  equally  weighted  observa¬ 
tions. 


Similarly,  if  m  unknowns  are  (indirectly)  observed 


h  times,  the  best  estimator  for  is  the  sample  standard 
deviation,  -5,  represented  by  the  expression  ("ref.  ll 


assuming  all  observations  are  equally  weighted.  The  value 
r  =  n-m  in  the  equation  is  known  as  the  "degrees  of 
freedom" . 
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A  set.  of  ft  observations  with  assigned  weights 

n 

u>,  j'-'j  is  equivalent  to  a  set  of  £  ^ i  observations 

lu 

which  are  all  equally  weighted.  Thus,  an  artificial  set  of 
observed  values  is  created  in  which  u)t  observations  are 
equal  to  the  1st  actual  observed  value,  observations 

are  equal  to  the  2nd  real  observation,  and  so  on  for  the 
remaining  weighted  observations. 

Given  an  arbitrary  set  of  weights,  the  set  may  be 
scaled  so  that  the  smallest  weight  has  a  value  of  ONE.  This 
scale  factor  is  known  as  the  variance  of  unit  weight,  . 

For  a  more  complete  discussion  of  this  topic,  see  Appendix  F. 

Therefore,  in  the  case  of  weighted  observations,  the 
best  estimate  for  the  value  of  the  standard  deviation  of 
unit  weight  is  given  by 


U-15) 


where  m  is  the  number  of  unknowns. 

In  matrix  notation,  eq.  (1-15)  is  written  as 


(I- 16) 
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the  number  of  unit  weight  observations,  is  given 


where  ft  , 
by  the  trace  of  weight  matrix 

h  *  trace  (  V)  *  2,  co,: . 

u»i 

.  XL 

The  standard  deviation  of  the  L  observation  (with 
weight  is  given  by 


G.  PRECISION  OF  ADJUSTED  VALUES 

The  elements  of  vector  X  (*i  ••  •  )  given  by 

X  *  (4V/0"  (M) 

represent  the  adjusted  values  of  the  unknowns.  The  matrix 
I  f  \-1 

14  ]  is  known  as  the  variance-covariance  matrix  Q,  and 

individual  elements  are  identified  by  the  term  . 

The  standard  deviation  of  each  adjusted  value  X£  is 

given  by 

(I- 19) 

where  j*  L  f  so  that  the  terms  are  diagonal  elements  of 

the  matrix  (A  W  A  ) 
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For  a  more  complete  discussion  of  precision  of  adjusted 
values,  see  Appendix  G. 

H.  THE  ERROR  ELLIPSE 

Position  errors  are  two  dimensional  and  must  be  evaluated 
in  terms  of  the  errors  along  the  x  and  y  axes.  Since 
the  maximum  and  minimum  errors  do  not  necessarily  occur  along 
these  axes,  the  orientation  of  these  maximum  and  minimum 
errors  must  also  be  considered.  Positioning  errors  may  be 
evaluated  in  terms  of  the  error  ellipse  (See  Fig  I-l) . 


The  greater  errors  occur  along  a  line  making  an  angle 
with  the  x-axis  (measured  anticlockwise)  such  that 


£.ot  1#B  ,  2^  ~  (1-13) 

Hu 


The  respective  standard  deviation  is  given  by  the  semi¬ 
major  axis  length  of  the  error  ellipse 


1  HtX _ 

~  D 


(  X- 24? 


wh^re 

D»  ■*  <•  1* 

The  smaller  errors  occur  along  a  line  perpendicular  to  5a  < 
and  the  respective  standard  deviation  is  given  by  the  semi¬ 
minor  axis  length  of  the  error  ellipse 


g  «*«.  Ht  Z 

•=)<• '  lit  t  p 


(  I  -  25  ) 


The  derivation  of  these  equations  is  presented  in  Appendix  H. 
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APPLICATION  OF  LEAST  SQUARES  ADJUSTMENT 


The  determination  of  a  vessel's  position  at  sea  is  a 
typical  hydrographic  problem  for  which  the  least  squares 
adjustment  is  particularly  well  adapted.  Various  methods 
can  be  used  for  fix  determinations.  In  this  thesis,  the 
following  three  methods  will  be  presented: 

a)  fix  by  N  azimuth  angles 

b)  fix  by  N  sextant  angles 

c)  fix  by  two  range  distances  and  one  azimuth  angle. 

For  each  method,  the  least  squares  adjustment  is 

applied  in  the  following  manner: 

1st:  solution  of  the  problem  for  particular  conditions 

2nd:  numerical  example 

3rd:  solution  of  the  problem  for  general  conditions 

4th:  formulation  of  an  algorithm  for  the  general 
conditions  case 

5th:  implementation  of  the  algorithm  in  Fortran  language 


( Nlortljii»9  ) 


Lu  c  C  5 


v 


N 


S 


(#■  *) 


X  (Easting) 


FIG  II-l:  FIX  BY  3  AZIMUTHS 
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FIX  DETERMINATION  BY  AZIMUTHS 
1.  Solution  for  Azimuths  from  3  stations 


a.  Determination  of  Adjusted  Coordinates 

Given  a  positioning  problem  as  diagramed  in 
FIG  II-l,  where: 

A  |p  —  is  the  observed  azimuth  from  station  1  to 
vessel  position  P 

A jp  -  is  the  observed  azimuth  from  station  2  to 
vessel  position  P 

AdP  -  is  the  observed  azimuth  from  station  3  to 
vessel  position  P 


)  -  are  the  grid  '-oordinates  of  station  1 

*  are  the  grid  coordinates  of  station  2 

(**,)*)  -  are  the  grid  coordinates  of  station  3 

the  grid  coordinates  (xy)  of  vessel  position  P 
will  be  determined. 

Step  1)  Formulation  of  observation  equations 

The  analytical  expression  for  the  azimuth  from  station  l 
)  to  P  (*,y)  is  given  by 


A*  ip  (ra4;3f»*K  -  £(*2). 

The  function  F(xy)  must  be  expressed  in  a  Taylor's 
series  around  an  "initial  position",  Pq  ,  whose  coordinates 
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are  defined  as  xg  and  y0  .  Evaluating  the  zero  and  first 
order  terms  of  the  series,  the  following  expression  is 
obtained : 

As  -  tan'1  *-*■•  +  >  ax  -  t  *>-*4  )  Ay 

ip  (xo-xl)1*  (yo-yi)4 


Designating  the  distance  and  azimuth  from  station  L  ) 

to  the  "initial  point"  (.*o, Y®)  by  SLO  and  A^t  io 
respectively,  then 

5 i.o  *  [  (  *0  -  xt  )*  ♦  ( y„  -  yi  )l  ]  /2 


and 

-  tan 


*D  -  *«: 
ya  -  Vc  • 


Therefore , 


y«>  -  xi 
‘Sio)1 


Ay 


and  the  observation  equations  will  be  (for  i=l,2,3) 


s 


$ 

l& 


B 


M 


where  ^  is  the  observed  azimuth.  In  the  matrix  form 
/4X-L-V  ,  the  obs.  equations  will  be 


(5 1*)1 

<S*0l 

(Si#)* 


Xo-x» 

(S,o)* 

(Sj#)1 

*o-Xj 

(5  30) 1 


ax 

dV 


3/»  - 


-4*10 

Ja 

/W»» 


r  i 

V, 


Vt 


The  angles  must  be  expressed  in  radians. 

Step  2)  Normal  equations 

Forming  the  normal  equations,  the  adjusted  values  for  fax 
and  will  be  given  by 

X‘(Ar'WA)'‘  (/Twl) 

ft 

where  ^  V 

4  V 

Step  3)  With  the  values  A  X  and  £y  a  new  "initial  point" 

Pi  will  be  obtained; 

x'0  -  to  +  A  X 
3 'or  +  A? 

and  the  procedure  may  be  repeated  in  an  iterative  way  until 
the  increments  &X  and  Ay  become  vanishingly  small.  Then, 
the  most  probable  values  for  the  coordinates  x  and  y  will 
coincide  with  the  coordinates  of  the  last  "initial  point"  obtained, 
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b.  Precision  of  Observations  and  Adjusted  Values 


i 


R 


I 

► 

l 


Step  1)  From  observation  equations  A  s  V  ,  where  JXT 
has  been  obtained  by  the  least  squares  method ,  the  residuals 
are  obtained,  i.e.,  the  differences  between  the  "true"  and 
observed  values  of  the  parameters. 

Then,  the  standard  deviation  of  unit  weight  is  given 
by 


•o  * 


%  2  2 
»,  V,  4  lu*  Va  4  uj5 

u ,  *  +  uJj  -  m 


Yl 


where  rr>  is  the  number  of  unknowns  observed.  For  that  problem, 
the  unknowns  observed  (indirectly)  are  x  and  y  (m*2) . 

Therefore,  in  matrix  notation,  the  above  equation  is 
expressed  as 


s0  ;  /  vTw  v 

v  trace  (  W) 


where  trace  (W)  -  uj,  -t  *  0)3  # 


Step  2)  The  standard  deviation  of  each  observation  (with 
weight  )  is  given  by 


li'Sti*)  . 
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Step  3)  The  standard  deviations  of  adjusted  values  are  given 
by 


The  covariance  is  given  by 

<  <  * 

a  IZ  •  '“'a  . 


Step  4)  The  correlation  coefficient  p  between  x  and  y 
is  given  by 


■5x3 
S*  Sy 


c.  Error  Ellipse 
Given  the  matrix 

Q,(Atwa)'‘.- 


91% 
9a  x 


and  the  standard  deviation  S0  of  a  unit  weight  observation 
the  error  ellipse  parameters  will  be  determined. 

Step  1)  Obtaining  the  value  D, 

D*  4  ]/z 


Step  2)  Semi-major  axis, 


2  in  in 
i<i 4  ii*  ■  D 


w 


Step  3)  Semi-minor  axis, 


S 


b  = 


^XX 


*=!«<  *  iiz  +  D 


A 


Step  4)  Determining  the  angle  (measured  anticlockwise) 

from  x-axis  to  semi-major  axis  : 

4.1)  The  angle  Va  will  satisfy 


tan  (2j0)  -  fan  XL  - 

^fii  ~  9 ii  • 


4.2)  For  computer  applications,  XL  is  defined  to 

fall  within  the  following  limits:  —  TT/2  *  <  tr/2 

Then, 

a)  if  ,  choose  y  =  TP/ 4 

b)  if  ^  and  j(L  ^  a  ,  choose  y-  Xl/i 

c)  if  #  in  and  _n_  ^  o  ,  choose  S'*  (-fl  +  TT)/^ 

4.3)  The  intersection  of  error  ellipse, 

1  t  c,  '2- 

^az  *  ”  2  2  x  2f  +  ^  "  y  "  CT*'  ^*1  -  ° 

with  the  straight  line  y  •*  xt»«y,  is  given  by  x,  and 


such  that 


^  (I  y-ZZ  So 
-  *  9ia  2T  +  <?..  W*  2T 


3? ,  x,1  .  fcd^> 


4.4)  Considering 


n a  *  1 

P.  »  x,  ■*■ 

D„l  ,  [IS*  t  Sk)  /i  ]* 


a)  if  D*  >  0*  ,  then  %0  ,  y 

b)  if  pj1  <■  p*  ,  then  %0  2  y  ti  /  2  . 

Numerical  Example 

a.  Determination  of  Adjusted  Coordinates 

The  U.T.M.  grid  coordinates  of  shore  stations 


FIG  ll-l,  are: 


COORDINATES  1  LUCES  (#1) 


MUSSEL  (#2)  I  MB 4  (#3) 


x  (EASTING) 


595,794.5 


597,967.8 


y  (NORTHING)  4,055,042.7  4,053,453.2 


603,425.2 

4,053,917.2 


For  illustrative  purposes,  the  standard  errors 
for  azimuth  observations  made  at  each  station  were  assigned 
the  following  values:  tf”,  r  o.  o  z  >  4j  s  o8.  » :  4  ,  4%  ,  0®  0  « 


The  observed  angles  at  each  station  were: 

P  -  LUCES  -  MUSSEL  =  =  50. *  164, 

P  -  MUSSEL  -  LUCES  *  <*x  =  99.*  360, 

P  -  MB 4  -  MUSSEL  =  <*3  =  47  *  865  . 

Step  1)  From  the  grid  coordinates,  the  following  azimuths 
between  stations  are  obtained: 

A,a  ,  4ai-T  £  l**-  /  C jt-iO  ]  *  <16?  (8 1 , 

A X |  a  A  +  16®  *  3® 4?  1 8  I  , 

»  tarC*  £(  *%)/  s  265?  14©  t 

Step  2)  The  observed  azimuths  will  be 

A\p  =  AIZ  -  o £  ,  s  76?  0(7, 

Aap*  Aa.t  t  o(z  ,  4*5?  S4I, 

Asp  «  A33.  f-  0(3  » 

Step  3)  Formulation  of  observation  equations 

3.1)  The  first  "initial  point"  is  determined  by  the 
intersection  of  azimuth  lines  from  stations  #1  and  #2  expressed 
by  the  equations 

f 

4  3  *  f*«  ( x-  x.) 

„  *  mx  c  x-**;  . 

Solving  these  equations  simultaneously  to  find  x  and  y  , 
these  values  are  used  as  the  coordinates  {Xo,yo)  of  "initial 


point",  where 


(X0  .  6  oo  }  877.  5 

\  ye  -  ^036J3oS.k. 

3.2)  Determining  azimuths  between  stations  and  "initial 
point"  fjj 

All0  -  tan"1  [Uo  -  X.J  /(yfl  -  y.  )]  s  76°  <5  17, 

42„  ,W [or.-vt)/<y, 45'j4t, 

^23a  .  tan  [u.  -X3)/  (Vo  -y»;  ]  =213*  ISS. 

Then,  evaluating  the  elements  of  the  L  matrix: 


*ie  - 

-  0.  odo 

-  0.0 

ra<i, 

^2ff  -  ^*1© 

-  -o*.ooi 

-  -  (5.  ooo 

ol7  5 

-  ^*3© 

-  -  CL  (80 

S  -0.003 

[4i  6 

ra<i. 

3.3)  Determining  squared  distances  between  stations 
and  "initial  point"  f©  , 

lo  )  s  (x©-x»)  +  ( y©  -  y» ;  »  27j438/895, 

LSl0  )2  ~  (v©  -  x*)1*  .  i6,*/8,  5ar, 

(S3£>  »  (*.  -y3)1*  Cy©-ys)*»  u,  zoe,6(3. 

3.4)  Then,  evaluating  the  elements  of  the  A  matrix, 

(y©  ->»  V-5 1©  -  o.ooo  1  -(v©->0/3»d  a-o.oooiesz 

tyo-Y*)/«$i©  ,  o.  ©ool7l9  -(y©-^2)/5ao  if -<a.o<3<5  I  75  ) 
(Y©-Y*)/  Szo  *°>  d©o  195*1  -(*©-*s)/S3*  a  0.006  2*61  . 


3.5)  Therefore,  in  matrix  form,  the  observation 


equations  are  written  as 


IS 

“ 

•  «■ 

r*  * 

0.000  066  1 

*  0.000185  1 

lAxl 

0.000  O0O 

v, 

0.000  171  8 

-  0.000  1751 

Wj 

-0.000  017  5 

X 

V» 

0.0  00  1454 

0.000108  7 

-0,005  1416 

V* 

- 

1 

- 

Step  4)  Solution  of  normal  equations 


4.1)  Determining  the  weight  matrix. 


<r,  * 

O.O40 

— »  1  /  (T7X  a 

2  5  00 

(fi  % 

0.024 

—  *  /85*  . 

1736 

«« 

t? 

0.  0  1  9 

— *  l  /  , 

3086 

Considering  the  least  weight  equal  to  ONE,  it  will  be  obtained 


that 

u)|  -  1.  A  4 

u)%  .  1.  oo 
u)3  -  1.7  8 


or 


144  o  o 
o  l-o  o 
o  o  ].78 


4.2)  The  solution  of  normal  equations  is  given  by 

X  *  (atw  a  f 1  LAt  W  l  ); 


4.2.1)  obtaining  matrix 


atw 


0. 0 Oo  066  4 
-0.000  166  7 


4.2.2)  obtaining  matrix 


0,000  )7/9  0.000  340  7  1 

*0.000  {75  /  a 00037/  5  J  , 


AtW/1  , 


0.000  000  10 
A 000  000  03 


0.000  000  05 

0. 000  000  1/5 


4.2.3)  obtaining  matrix  Q  =  (/i'W  A  ) 


•  t 


Q. 


to,  A7»,  2  o<* 
i;  510^  68 1 


<;S7 o,  68 1 

5,  235,  602 


4.2.4)  obtaining  matrix 


aV  l 


0.  000  0  0(  I 
0.000  001  2 


4.2.5)  finally,  vector  X  is  evaluated  by  solving 
the  normal  equations  : 


X. 


-  8.6 

-4.6 


Step  5)  First  adjusted  values  of  x  and  y 

With  the  increments  A*  and  Ay  a  new  "initial 
point"  is  obtained: 

f  Xo  =  6oo,  877.5  -  3.6  s  6oo,  867.8 
1  yoz'ijOS 6,3<>8.4-  4.6  ,4,056,303.8  . 

Step  6)  With  the  new  values,  for  the  "initial  point",  the 
procedure  indicated  in  steps  3.2,  3.3,  3.4,  3.5,  4.2  and  5 
is  repeated,  and  with  the  values  now  computed  for  AX  and 
An  a  "closer"  initial  point  is  obtained. 


Step  7)  This  procedure  must  be  repeated,  in  an  iterative  way 
until  the  increments  A*  and  become  vanishingly  small, 

or,  in  practical  terms,  converging  to  within  a  specified 
tolerance.  Then,  the  last  "initial  point"  obtained  will 
coincide  with  the  most  probable  position  for  P(x,y). 

b.  Precision  of  Observations  and  Adjusted  Values 

Step  1)  The  residuals  are  obtained  introducing  a  x  »  -3.6 
and  -4.6  into  the  observation  equations.  Therefore, 


0,0 oo  4 ©6  4 
-o.  ooo  9  26  3 
O.ooo  3 oo  3 


Step  2)  Obtaining  scalar  y  W  V , 


VTW  V  *  o.  ooo  oo\  oSS  . 


Therefore,  SQ  =  0.000  6992  radians  . 


Step  3)  Obtaining  standard  deviation  of  each  observation  , 


5.* 


Then, 


S  ,  .  o.  ooo  7  ra<i 
5 2  -  O.ooo  6 33 Z 
Sj  .  o.  ooo  ra<* 


*  0.° 03  3 

s  o* oq  o 
s  0*03  O  , 


Step  4)  Standard  deviation  and  covariance  of  adjusted  values 


j 

I 

I 

! 

I 


x  and  y  , 

Sx  =  5  2.26 

$7  =  5  i  1 .60 

Sjt y  s  ^it  s  -0. 768  , 

Note,  Sy  and  Sy  are  expressed  in  the  same  units 
as  the  grid  coordinates . 

Step  5)  Correlation  coefficient, 

P-  -  _  0.21  # 

'  Sx  • 

c.  Error  Ellipse 

Given  : 

So  =  o.  ooo  6  99  Z 

sr  1°;  ^»7(j  2o£f  , 

=  S,  6o  2. 

^ia  *  -  f  * 

the  error  ellipse  parameters  will  be  obtained. 

Step  1)  Determining  D, 

D--  [(3,.  -  ]  .  6,  «■*,  7°3  • 


Step  2)  Semi-major  axis  , 


Sa=  50 


2 ill 


D 


'A 

C  1.36  . 


Step  3)  Semi -minor  axis  , 


£l  5  S, 


Z  <3n. 

%•  +^ta  *  0 


*  1-57  . 


Note,  Sq  and  sb  are  expressed  in  the  same  units 
as  the  grid  coordinates. 


Step  4)  Determining  the  angle  q0  (measured  anticlockwise) 
between  x-axis  and  semi-major  axis  5 3 

4.1)  The  solution  of  equation 


tan  JCL  -  1  - 

^  «4  - 

is  Jl  «  -  3 of  56  4  . 


4.2)  Since  cj,,  andjl<0  ,then 

VA  7<t°52  . 

0  2 

i  2 

4.3)  Obtaining  X,  and  , 

.  *i%  0.  I7S 

y,!s  2.«8 . 

4.4)  Obtaining  and  Da*f 

(  Pi  s  y 1  ♦  3‘  -  2.  4  6 

t  P?  -  *  SUJ  /Z  j\-  3.3*  . 

4.5)  Since  ,  then 

g'a*  #*5$%  >6^52  (See  FIG  II-2)  . 
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FIG  II-2 :  ERROR  ELLIPSE 


All  of  these  computations  may  be  compared  with 
those  shown  in  the  computer  output  section  on  page  148  . 
Differences  in  the  results  are  due  to  the  fact  that  the  cal¬ 
culations  illustrated  on  the  preceeding  pages  were  only  carried 
out  for  one  iteration. 

3 .  Solution  for  the  General  Case 

The  solution  will  be  presented  in  a  way  that  can 
easily  be  implemented  by  an  algorithm  satisfying  a  modular 
design. 


a.  Determination  of  Adjusted  Coordinates 


1 

►  '» 

!' 

I 

1 

.** 

*  * 

*  ' 

3 

i'l 

L\ 


% 


H 


Given: 

a)  the  grid  coordinates  of  N  statons  Si  L*isc  ), 
where  the  x-coordinate  represents  EASTING'S 
and  the  y-coordinate  represents  NORTHING'S 

b)  the  azimuths  Aw  p  N)  from  stations 

5^  to  vessel's  position  P  (  ) 

c)  and  the  standard  deviations  OV  (  i  *  i,*  ) 

of  observed  azimuths, 

the  adjusted  coordinates  for  P(x_y)  will  be  determined. 
Step  1)  Weight  matrix  w 

1.1)  Squaring  the  inverse  of  standard  deviations  , 

s  -  »  (  &  »  tj  1 , |  N  )  , 

1.2)  Designating  by  u>'if  the  least  ,  the 

weights  u^;  will  be  obtained; 

C  L  *  1  i  *■)•••  j  **  )  , 

uj'k 


that 


1.3)  The  elements  of  square  matrix  W  will  be  such 


0  ,  P»r  C  ^  j 

***«:  >  for  L  -  J 


( c)  J 


',2, 


N  ) 


Step  2)  Observation  equations 

2.1)  Determination  of  first  "initial  point" 

2.1.1)  Designate  by  A\(p  an  observed  azimuth  Ai,p 
(i=2,3, . . . ,N)  such  that 

tan  . 


If  no  such  azimuth  is  available,  then  the  vessel's 
position  is  undetermined . 

2.1.2)  The  intersection  of  the  azimuth  line  Afcp 
from  station  K  )  with  the  azimuth  line  Aip  from 

station  1  (  )  determines  the  first  "initial  point" 

P0  (Xo,7©  )  •  Therefore, 

a)  if  (n=0,l)  ,  then  will  be  given 


by 


X0  - 

(  22*  '  &K?  )  •  *K  )  y 


b)  if  A(cp  *  h  TT  (n  *o,  t) ,  then  will  be  given 

j  X*  *  y  (r 

y9*  yx+  -  ^ip  J  *  (*«■  -  *»)  , 


c)  otherwise, 


Po 


will  be  given  by 


1  X0  «  -t*  "  *  rWl  ~  m*  ** 

J  »"i-  ^ “ 

[  zt x*  rn*.  ^  xo  -  x j ) 
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where 

mis  tah[(5tr/n-  Aie  ] 
me.  tan  -  A  kb  ]. 

2.2)  Determination  of  azimuths  A*it  between  stations 
and  "initial  point"  Pa 

2.2.1)  Two  angles,  At  and  satisfy  the 

equation 

Aii0'ia»'  U*> 

^So  "  3  i 

Also,  Ai  IQ  must  be  a  positive  angle  between  0 
and  2 T.  Since,  in  general,  calculators  give  a  solution 
between  (**  TT/2)  and  (+TT/2  )  ,  a  criterion  will  be  estab¬ 
lished  for  selecting  the  valid  solution. 

2.2.2)  Criterion  : 


a) 

if 

>  *  3  * 

and 

y<>>  t  then 

to*  ir/z 

b) 

if 

and 

f  then 

=3"/' 

c) 

if 

and 

2©  >**.'  /  then 

y^Zto  *  o 

d) 

if 

XosXL 

and 

<?c  ,  then 

A*io  =  t7 

For  X©  XC  and  ^  ,  designate  by 

the  solution,  given  by  a  calculator,  of 


oC;  - 


-  ta  h“'  m  *  i 


(  i  s|,2 .  A/  )  . 


Therefore, 


e) 

if 

<■»>« 

and 

X©>x*; 

,  then 

~  ^ic 

f) 

if 

°<c»<® 

and 

*©> 

,  then 

A'iio 

g) 

if 

O^c.  >0 

and 

Vo  <XC 

,  then 

A  *£o 

h) 

if 

and 

v*  <  yt- 

,  then 

An*  =^2ir 

2.3)  Determination  of  elements  of  matrix  L* ; 

U*  Ait- 

2.4)  Determination  of  squared  distances  between  $  i  (*£  $i) 
and  £(*•»>): 

1  't  t  t 

Sio  a  +  (  L  *  t,2  A/  )  , 

2.5)  Determination  of  elements  a^j  (l*  Ij"-j  Af  '  J»  f,3  1 
of  matrix  A  • 

at  .  y° ~ 2L  u  .t.z . w  > 

(  Sc’.)* 

a{1.  .  x»-*< 

Step  3)  Normal  equations 

3.1)  Determine  matrix  A  W  (a  matrix  2  x  N), 

3.2)  Determine  matrix  ArWA  (a  matrix  2x2), 

3.3)  Determine  matrix  l (a  matrix  2x2). 


* 

Since  is  a  symmetric  matrix,  then  its 

inverse  matrix  will  be  Q  =  ,  also  symmetric, 

such  that 

^  *i  r  r 

A.  A,t  Q„  Q,»  -  1  o 

,  Aai  .Qji  Qit,  ^ 

It  can  be  shown  that 

Qii  a  —  1  (  ”*  .  A\%  ) 

t?n  «  On  «  A  14  /  [At*  -  A ii  . 

G^it  *  — 4  «♦  /  <4%% )  . 

3.4)  Determine  matrix  ^rW(.  (a  matrix  2x1). 

3.5)  Finally,  determine 

-X’.  dx  ]  c  MTWi*r'  (Ar W  L). 

w 

Step  4)  First  adjusted  values 

With  the  values  and  the  coordinates  of  the 

new  "initial  point"  Po  (  x*o  ,  J  <»  )  are  obtained: 

* 

«  )(#  *  +  4  ^ 

.  . 

Step  5)  2nd  iteration 

For  obtaining  a  "closer”  initial  point  repeat  the 
steps  2.2,  2.3,  2.4,  2.5,  3,  and  4. 
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Step  6)  Next  iterations 

Repeat  Step  5  until  and  6^1  become  vanishingly 

small,  or,  in  practical  terms,  converging  to  within  a  speci¬ 
fied  tolerance. 

Then,  the  adjusted  values  for  x  and  y  will  coincide 
with  the  coordinates  of  the  last  "initial  point"  obtained, 
b.  Precision  of  Observations 

Given  N  (number  of  stations)  and  the  matrices  A,  X, 

W  and  L  determine : 

Step  1)  Matrix  of  resid”*is  V  (a  matrix  Nxl) 

v, A  X-  L 

Step  2)  standard  deviation  $0  of  the  unit  weight  observation 

2.1)  Obtain  VTWV  (a  scalar). 

2.2)  Obtain  trace  of  weight  matrix  W; 

V 

t race  (  W  )  »  ^ 

,  # 

u*i 

where  is  a  diagonal  element  of  weight  matrix  W . 

2.3)  Finally,  $0  (in  radians)  will  be  given  by 

s.*  /  . 

H  tra c«(W)-2  ’ 

Step  3)  Standard  deviation  of  each  observation  (with 

weight  ) , 

$  t  *  - — — ° . —  (.£•»  i  jt/.MjN) 

where  5,1  is  expressed  in  radians. 
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FIG  II-3 :  FIX  BY  3  SEXTANT  ANGLES 
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B.  FIX  DETERMINATION  BY  SEXTANT  ANGLES 


1.  Solution  for  3  Sextant  Angles  (Between  4  Stations) 
Given  a  positioning  problem  as  diagramed  in  Fig  II-3 
in  which: 

cLl£2,  “  t*ie  observed  sextant  ingle  from  P  between 
stations  1  and  2 

**2£3  “*  observed  sextant  angle  from  F  between 

stations  2  and  3 

0*^^  -  is  the  observed  sextant  angle  from  f*  between 
stations  3  and  4 

and 

)  **  are  grid  coordinates  of  station  1 

“  are  grid  coordinates  of  station  2 

'  -  are  ?rid  coordinates  of  station  3 

(X^A)  -  are  grid  coordinates  of  station  4 

the  grid  coordinates  (xy)  of  a  vessel's  position  P 
will  be  determined. 

Step  1)  Formulation  of  observation  equations 

The  analytical  expression  for  the  sextant  angle 
u  f  (ifl)  •  from  the  vessel's  position  P  (x,y)  ,  between 
stations  L  (*£,}£)  and  L.+ 1  (  ,  XV* )  is  9iven  by 

^LpCu+O  ^,n  lra^,ans)  5  I)  “  ~ 

tan*1  ZiLLzli.  -  tan 1  Jii  *  iL,  *  Fc*>y). 
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The  function  F(x,y)  must  be  expressed  in  a  Taylor's 
series  around  an  "initial  position",  ,  whose  coordinates 

are  defined  as  and  y0  .  Evaluating  the  zero  and  first 

order  terms  of  the  series,  the  following  expression  is 
obtained : 


Rx3  )  *  Ai  e  -  Aifi  a 


tan*'  Y»'  ~  *«  .  tan'  + 

3  i-n  -  21  a  2Jl  ~  2!  o 


-  3  i+, 


3#- 


i  2f*  ~  21  i+i)*  +  m*L+*  (y#  *  yi  Cx*  -xi)1 


4* 


-  Xc _ 

cy^-y^1  *  (x.  -Xc  >* 


Xq  *  Xtf« 

cyo-ytv.)**  i*o 


Designating  by  5*£  and  S^^Vo  •  and  ^oi  and 
the  distances  and  azimuths  between  "initial  point"  (Xo,2f») 
and  stations  L(Xt,20  and  ^  )  '  respectively, 

then 


..  r  *  1 1  '4 

5^1  s  L c Xt  —  x»)  •♦■Cvt-  J 

^o(i*  i)  s  [*C  "  ** )  +  C  ytVi  *  1 

/4?0(;  .  tan-4  [(x; -xo)  /  -  >)  ] 

s  ^n~*  £  (Xu*' ~  ^  /  c  y»;+«  “  y* )  1 
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and, 


if  (i*0  •  * 

[(3.-3^.  )/(5fl  )*-  (3.-  )/  £  5 o  £  >*  ]  4X  + 

[lx. -XI) /(Sot  )l-  (X,.  . 


Therefore,  the  observation  equations  may  be  expressed  as 
(for  i  =  1,2,3) 


V;  * 


^  l<»» 


-  X* 


(Sou..))1  (SoO*  J 

» 

o I. 


+ 


Wo  -  Xt.  •  Xu* I  ^ 

.  (SoO1 '  csa(t  «,))*] 


i.eiin)4  A*o(c< 


where  ( t+0  ttle  observed  sextant  angle.  In  matrix 

form,  the  above  equation  is  expressed  as 


V* 

where 


A 


/IX- 

L 

X<,-Xi 

<kSoi  )* 

CSo,)% 

y»  -x*. 

C5oi)‘ 

*e-1l 

Xo-^S 

(Saa)1 

(i.»X 

Vo  »  Yt 

(Sa>)’ 


($«*)’ 


»xt 

CSai)* 

X»  ~  *3 
(So3)Z 

x>  **x> 
rSoO’ 
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rA 


V- 


In  that  result  it  should  be  noted  that  the  angles 
are  expressed  in  radians. 

Step  2)  Normal  equations 

Forming  the  normal  equations,  the  adjusted  values  for 
AX  and  Ay  will  be  given  by 

X  *  (AJ WAf‘  u’wu. 

Step  3)  With  the  values  Ax  and  a  new  "initial  point" 

P  o  (*0,^0)  is  obtained, 


(  *  a  *  *o  + 
1  + 


and  the  procedure  will  be  repeated  in  an  iterative  way  until 
the  increments  AX  and  Ay  become  vanishingly  small. 

Then,  the  most  probable  values  for  the  coordinates 
x  and  y  will  coincide  with  the  coordinates  of  the  last 
"initial  point"  obtained. 

2.  Numerical  Example 

Referring  to  FIG  II-3,  the  U.T.M.  grid  coordinates 


I 


of  shore  station  are  : 


Coordinates 

X (EASTING) 
y (NORTHING) 


MB 4  (#i) 
600,425.2 
4,053,917.21 


USE  (#2) 
600,372.0 
4,051,216.9 


MUSSEL  (#3) 
597,967.8 
4,053,453.2 


LUCES  (#4) 
595,794,5 
4,055,042,7 


The  observed  sextant  angles  are  equally  precise;  thus,  the 
weights  will  be 

CO ,  i  c 0-2,  s  0)3  :  1  , 

and  the  weight  matrix  W  is  the  identity  matrix 


o 

o 

l 


The  following  sextant  angles  were  measured; 
MB 4  -  P  -  USE  =  of  jp  j,  =  49*.  927 

USE'-  P  -  MUSSEL  =  0<2f3  =  38*.l30 
MUSSEL  -  P  -  LUCES  =<*3,4  =  30#.  396 


Step  1)  Formulation  of  observation  equations 

1.1)  Determination  of  first  "initial  point"  Pa  (Xo,3o) 

1.1.1)  The  first  "initial  point"  will  be  the 

point  determined  by  the  two  sextant  angles  and 

such  that 


I  P  2 


X7  —  *© 


oL 


X© 

3,-  J* 


and 


Xs-xo  _  x*-x0 

^i>yg _ y» -*« 

j*  Xj  -  x#  y%  ~y» 

*»-*>  . 


1.1.2)  After  some  algebraic  manipulation,  it  will  be 
obtained  that 


■j  e  a  C  *0  *  O 

where 

*  x,  -  X3 

C  -  t**  »<if  i  tand»E3 

*  -  *>-  -  *■-*.  -  3,  ♦  3* 

and 

3,  *» 2  *'  -  **-**"  illS  +  1*  •«  X*  *3  -  X;  x7 

JX  tan  <*,<,?  i»"  •'l** _ _ _ 

**  -  **  -  X<  ~  X*  -  -  J  I  +  jf  3= 

ta»  o^ej  tdr,  oi  tp2 

1.1.3)  The  value  X0  will  be  a  solution  of  equation 

UK,!t  8/,  *5:0  (I  -i) 
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where  U  mta««,rl.  ic4*  i  j 

R  =.  tan  t  [iCD- dyt+yi)- (Xl*  xz)1^ - 

Cfri-Vi)  +  3,  - 

and  r  i  7 

5  »  ia«o(|f2  [  D  -  D  ( ^ I ♦  3a)  t  X\  *x  +  y,  3a  J  - 

&[*,“  *»)-?/  *a  4  ji  x, 

1.1.4)  IWo  solution  sets,  (X#l  ,;}«,)  and  (x6t,ja2), 
satisfy  eq.  ( II— 1 ) .  The  valid  solution  corresponds  to  the 
solution  set  that,  introduced  into  the  following  expression, 
yields  the  value  that  best  approaches  tan  0(^2.  • 


(yz-y*) 

— - - - "‘Van 

(?!-?»)  (3,  -».J  +  (*i-  X.)  (y,-  yi) 


=<  C»-z) 

^IP  z 


1.1.5)  Using  the  numerical  values 

tan  oCj  =  1.1887  tan  ^^,£3  ~  0.7850 
X, =  603,425.2  ^,=4,053,917.2 

Xz=  600,372.0  y2=  4,051,216.9 

y^=  597,967.8  4,053,453.2 

it  will  be  obtained 

C  =  il. 109096  D  =  -2,613,387.9 

U  =  147.885  403  R  =  -1.776434  x  108 

S  =  5.334  734  3  x  1013  . 
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The  two  solution  sets  satisfying  eq  (H-l)  are 

=  600,833  and 

1  =  4,056,325 

f  *oz  =  600,390 

l  ^1  *  4,051,405 
Introducing  the  first  solution  set  (y0l  )  into  expression 
(II-2)  the  value  1,2923  will  be  obtained.  Introducing  (*Vi, 
y^x)  into  the  same  expression , the  value  -0.9944  is  obtained, 
Since  the  first  set  is  the  one  that  best  approaches  the  value 

of  ■ban  oL  -  1.1887,  the  coordinates  of  first  "initial  point' 

If* 

are 

I  *  600,833 

l  ^  =  >1  =  4,056,325  . 

1.2)  Determination  of  azimuths  between  "initial  point' 
and  stations  Sc  ( *£,2^  /  (i=l,2,3,4)  j 


Alo  1  =  tan'  =  >32.°«ea 

Aloz  =  tan*'  [lx,-*;  /[  x  - y.)  ]  =  lesh si 

A?oi  ,  tan”'  [u3-*.)/  ]-  2Z4f634 

Alal,  r  t an1  )/{yu-  >)  ] -255?  721  . 


Then, 


o<tez 

* 

-  Ai  at 

-  1*7 HZ 

=  -  0.  O<40  0JF  t 

rad 

* 

fit  at 

i  -/*  iki 

?  -  0. $10  705  5 

had 

*30* 

t 

fit 

-  Aiok 

*  •  0.  35 1 

r  -  0.  Oof>  0  20  l 

1.3)  Determination  of  squared  distances  between 
and  stations  ; 


lS0l  )*- 

+ 

2 

*  12,517,002 

[  S  ox  1  =  ( ) 

+ 

=  26,305,207 

(  6 03  1  = 

+ 

(y.-n.)1 

-  16,456,606 

( 5^a  )  = 

+ 

i#.  )* 

=  27,030,776  . 

1.4)  then , 


CS*t)' 

*jf  *^i  X  0-  OOO  00f  8 

(  So#)* 

Xo-xi 

—  0.000  22k 

<W 

^•"3*  7  -a.oooa/*  7 

csj' 

Xo-Xi  _ 

•x.-^  „ 

—  O.  D0O  /  A 

.  _  o.ooo  Uf  l 

Xo-x'v 

—  0.  0O0  0/2  3 

is**F ' 

IS")1 

cw* 

1.5)  Therefore,  in  matrix  form,  the  observation 
equations  will  be  expressed  as 


4  0.  000  00/  8 

-  0.000  2  24  6 

*  — 

-  0.0*0  07 5  6 

r  - 

vx 

-  o.  00O  0l«f  7 

-  o,  ooo  zf6  3 

4  * 

4  * 

-  028  7  MS  6 

- 

Vt 

-  0.  000  117  J 

-  o.  ooo  o  / 1  3 

u  ** 

—  0 . 006  82*  2 

to  j 

.  *  to 
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Step  2)  Normal  equations 


The  solution  of  normal  equations  is  given  by 

X  =  lAr  W  A )"'  MTWL)  i 

2.1)  obtaining  matrix  AT  W, 

A  W  S  0-000  OOI  0  -O.OOOO101  -0.000  It 71 

-d.Ooo  12k  6  -o.oootS66  -o.ooooiid 

2.2)  obtaining  matrix  atwa  , 


T  —  8  — 

AWA.  i  .6  s  >*  i  o 

y  1C9  7.  5-1  y  i<?‘ 

2.3)  obtaining  matrix  Q  =  [4TW4)  *, 


-  ■ 

Qa  6«,<«*8,178  -3,^72,073 

-^<i7i,  13, sn, 606  # 

2.4)  obtaining  matrix  A  W  L  , 

ATWL  s  1.3*0  y  id"6 

.  (.371  ,  IO-5  J  , 

2.5)  finally,  vector  X  it  will  be  obtained 

Xz  35.8 

.131.3  .  . 


Step  3)  First  adjusted  values  of  x  and  y 

With  the  increments  AX  and  a  new  "initial 

point"  is  obtained; 

J  Xo=6do, 833  +  35.8  *  6oo, 868.8 

I  >  *  4,0$6;325  ♦  191-3  =  UjOS6}5o 6.  3  . 


Step  4)  With  the  new  values  for  the  "initial  point",  the 
procedure  indicated  in  steps  1.2,  1.3,  1.4,  1.5,  2  and  3  is 
repeated,  and  with  the  values  now  obtained  for  and  A  y 

a  "closer”  initial  point  is  obtained. 

Step  5)  That  procedure  must  be  repeated,  in  an  iterative  way, 
until  the  increments  AX  and  Ay  become  vanishingly  small, 
or,  in  practical  terms,  converging  to  within  a  specif ieu 
tolerance.  Then,  the  last  "initial  point"  obtained  will 
coincide  with  the  most  probable  position  for  P(xy). 

These  computations  may  be  compared  with  those  shown 
in  the  computer  output  section  on  page  149  .  Differences 
in  the  results  are  due  to  the  fact  that  the  calculations 
illustrated  on  the  preceeding  pages  were  only  carried  out 
for  one  iteration. 

3.  Solution  for  the  General  Case 

The  solution  will  be  presented  in  such  a  way  that 
easily  can  be  implemented  by  an  algorithm  satisfying  a  modular 
design. 

Given: 

a)  the  grid  coordinates  of  M=N+1  stations  Si  (  , 

ordered  in  a  clockwise  sense  around  vessel's 
position , 

b)  the  N  sextant  angles  d  L  e  Ci*l)  between  stations 

5  £  (  X£,  )  and  ^  C-t  i  (  ^C*i  '  ^  ^  t 
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a 

i 

■ 

3 

I 


c)  and  the  standard  deviations  0 c  (i=l,2,...,N) 
of  observed  sextant  angles, 

the  adjusted  coordinates  for  P(x,y)  will  be  determined. 

Step  1)  Weight  matrix  W 

Obtained  as  indicated  on  Step  1  of  subsection  II.A.3.a. 

< 

Step  2)  Observed  equations 

2.1)  Determination  of  first  "initial  point"  S?0  {X0,^a) 

2.1.1)  The  first  "initial  point"  (  *b,jf«)  will  be 
the  point  determined  by  the  two  sextant  angles  o/jpg  and 

such  that 


y-i  -  x0  Xi-  xo 

tan  <*IPJL  z  tan  (AlM-  Alol  )  »  — ^  - P'-ritL- 

(  +  Xx  -  Xo  #  Xt-  *o 
3a  “  So  3.-3* 


i 


and 

)^3  -  Xo  X2  ~  Xo 

^3  -  3z-  J* 

1  +  *»“  *•  .  x0 

?3  -  %  ?2  -  3*  . 


4 


2.1.2)  Test  for  undetermined  initial  position 
(See  FIG  II-4) . 


4 


?o  tXe>) 


FIG  II-4:  UNDETERI1I NED  FIX  BY  2  SEXTANT  ANGLES 

If  ,S,,  5Z  and  S3  belong  to  the  same 
circumference,  then 

y3=lSo  -  (  C*  ip*  ♦  ). 

Therefore,  the  dot  product  of  vectors  U  and  V 
(FIG  II-4)  will  be 


u.  v-  (/3-  *x)  f  (  -  |U|.  171  co&fi 


So,  the  condition  for  an  undetermined  fix  by  two  sextant 
angles  will  be 

-ft) _ 

(3l-ja,tJ[lvy/  +  fy.-s*')1]}  1/1  _ 

2.1.3)  Of  the  initial  position  is  not  undetermined, 
then  its  coordinates  can  be  obtained  as  follows: 

2. 1.3.1)  1st  case:  j  9$  and  o^2g3  £  9o* 

Let 

A  &  tan  oi  i  £  ^ 

0  »  tan  o*af  3 

F  a  A  1  4  B]  *  *•  “  *9 

6,  (:M*-**0//U  (/*  *-*3  30/3  ♦JfcX*  -I-** 

a)  If  F=0  ,  then 

G/F  r  P. 

Let 

Um  A 

R*-A(x,+  *»)  o.Oz 

5-  A[d- PC****)*  *iX*.  +3.^1  3+  PC*z-x,)-f  ,  . 

Then,  from 


LTy0% 


^  X©  -f-  _S> 


o 


obtain 


b)  If  E=0,  then 

V0  -  -  G/f  =  H  . 

Let 

v.  /A 

R  -  -  A  (^,  -*•  *«)  -  *1  +  x* 

5-  A  f^-Wlx,4  x2)+x,y»4  W(jr^)+  *'2z~ 

Then,  from 


obtain 


From  solution  sets  ■  >  'j* i  )  and  ^  choose  the 

one  that  best  satisfies  equation  (II-3) . 

c)  If  E  ^  0  and  F  +  0  ,  then 

:  {  F  /  E )  Xe  +  G  /  £  -  D. 

Let 

v=  A  (c% 0 

R  ,/l[2  CD-  Ciy,*3t)  -  (*,+  *»)]  -  C*  +  C  *2  -2x  + 

Sr  A[pl-P  (},♦**)+  *»>i  +  3i^l  +  -*»*»  . 

Then ^ from 

VjfcSR  x.  +  S  s  o  obtain 


z  u 

From  solution  sets  (x^,  ,  ytfl  )  and  (xfll ,  yol )  choose  the 
one  that  best  satisfies  equation  (II- 3) . 
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0  o 

2. 1,3. 2)  2nd  case:  °*\£2  *  90  and  ^2?!^  90 

Let 

B  X  tart  dZ£ 3 

E  «  $  c?»-  :>• »  +  *»-*» 

F*  B  -3a 

£>  *  B  O*'?*  -y*1%  +  **  :U  -  *»  dx  . 

a)  If  F  =  0  ,  then 

c  /  £  *  p  . 

Let 

R=*  —  ^  tf,  •+•  X») 

■S-PJ,-P(3i+3»)  +3,3,  +  *i  *»  . 

Then , from 

xJ"  4-  R  Xo  f  5  *  ©  obtain 

*ai  -Rt  >/«*-  *s 
2 

From  solution  sets  ^y0  )  and  choose 

the  one  that  best  satisfies 

“tan  <^ipi  :  (**■-*•?  (u-4) 

(djo.)0-  %>  +  (**-*•>  (**-  • 

b)  If  E  «  0  ,  then 

X0  ?  -  S/F  =  H  . 
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Let 

R*  -  ( 2\  + 

£>s  (41*  .  H  **)  ■+  *■  *x  t  3, 3*  . 

Then,  from 

+  R  +  S  »  O  obtain 

2. 

From  solution  sets  )  and  (y*,  ^*i)  choose  the  one 

that  best  satisfies  equation  (II-4)  . 

c)  If  E  ?  0  and  F  f  0  ,  then 

■a0, (F/ E)*.*  (  G/r>  «■  p. 

Let 

Us  C- H 

R*2  GP-C  (},*:>*)-  (*1  +  **) 

3  =  D*  -  P  U,-OiJ+  +  *»  **  . 

Then,  from 

2, 

Xf  X0  +  R  Ko  +  S  *•  obtain 

>/g*-4U3. 

2.  O 


From  solution  sets  Cx^,  )  and  (Xoirjat)  choose  the 
one  that  best  satisfies  equation  (II-4) . 


2. 1.3. 3)  3rd  case:  =3o  and  &  igz  ^  ^0 


A  s  tan  |g  2 
E  r  /4  0‘-»  +  *•-  ** 

F  *  A  ( X*  -*•)  +  *»  -  3*- 

^  *  4  (  X|  Xi  4  I  ^)x  -  yt*  x3  ’ 


3*  ?a  )  +  *.  -^1 


a)  If  F  =  0  ,  then 


^)o  »  & /  £  —  P  . 


R-  -  c *. ♦  *») 

5jP‘.  J!(  ?i  +  3».)  +  Si  3*  +  *i  x* 


Then,  from 


X>  4-  R  x*  +  S  -  o 


obtain 


v  -  Rt  \/R*-  4  S 

*0  S  i 


From  solution  sets  {Xoi  >2#)  and  (X0t  J«)  choose  the  one 
that  best  satisfies  equation  (II-3) . 

b)  If  E  »  0  ,  then 

Xo  r  -  6  /  f  -*  M. 

Let 

R  -  -  (  ^4^3/ 

5  r  W4-  H  (^4^)4  . 


Then ,  from 


Jo  +  &  Jo  +  S  s  o  obtain 

3o*  -Rt  v/g1- 

From  solution  sets  (X?  ,  'jty )  and  (*d,^oz)  choose  the  one 
that  best  satisfies  equation  (II-3) . 

c)  If  E  ^  0  and  F  ^  0  ,  then 

s  (F/e)y„  +  (  <=*/£)=■  cx«  ♦  d. 

Let 

U3  C  4.1 

Cp-C-C  jt+^3)  -  (  X*-* 

5r  P’-P  +  3*.  3?  . 

Then,  from 

VJ  Xo  +  P*  "Vo  +  3=0  obtain 

X..  -Rt\I Rx-  W5 

2  \J 

From  solution  sets  (>#I  )  and  (*©*,  , 2«i.)  choose  the  one 

that  best  satisfies  equation  (II-3) . 

2. 1.3. 4)  4th  case:  0^ =  90°  an<*  °^2F3~ 

Let 


£  *  Ui- 

F  =  *»-  *» 

6sxiyaOi^-^^-^  ^  . 
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a)  If  F  =  0  ,  then 

&/E  *  P  3nc i  Xi  , 

b)  If  E  *  0  ,  then 

X0  s  -  G/P  r  H  anA  y0  z  y  i  . 

c)  If  E  ^  0  and  F  ^  0  ,  then 

*  (  F/e) *o  *  (  G/E )  -  C  x„  +  P  . 

Let 

U*  C%  i 

Rs2CP-C(^4^)-  (.**.  +  *») 

6a  P*-  P  ^3  +  . 

Then,  from 


U*  X*  -J-  RXo  f  S=o  obtain 

xi,  .  -  Rt  V  R*-4US 
1  \S 

From  solution  sets  {Xoy  ,y0l  }  and  (•*<9t.f}dz)  choose  the  one 
that  best  satisfies 

(>-?•;  ( 3i-  ?*)  *  {Xz-X*)  C*l~  X°)  -  °  .  (31-5  ) 

2.2)  Determination  of  azimuths  A?o  L  between  "initial 
point"  (*«>}#)  and  stations  Si  (  *i  JC  ) 

2.2.1)  Two  angles,  and  /4l^  +  180  /  satisfy 
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the  equation 


=  ‘tan’1  (iai;2r,.  a/+«). 

yc-yc 

Also,  A* 0  ^  must  be  a  positive  angle  between  0  and  . 
Since,  in  general,  computers  give  a  solution  for  the  above 
equation  between  (-tr/2.)  and  (4?r/*r)  ,  then  a  criterion  will 
be  established  for  selecting  the  valid  solution. 

2.2.2)  Criterion  : 

a)  If  >  and  ,then 

b)  If  and  ,then  A^oi~  ^’/z* 

For  /  'ju  designate  by  oi0^  the 

solution  given  by  a  computer  of 

otoi  s  tah*  *  v,  ~ *».  Al+|)  , 

Vi-Jo 


Then : 

c)  If  O  and  XoJXc ,  then 

s  <*oL  +  • 

d)  If  o/0(,^o  and  X0  £  Xc  ♦  then 

A^Cti  S  °^oC  . 

e)  If  <0  and  Xo>  X  £  3  then 

a  +2.7T . 

f )  If  oIqL  <0  and  ,  then 

ql.  -  "*  TT , 
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2.3)  Determination  of  elements  LJ  of  matrix  L  : 

*<££(£««)  ♦  At  A  £  -  At0  £  t%|)  (  i  %  1^  ...J  Al)  , 

Note,  Ly  must  be  expressed  in  radians. 

2.4)  Determination  of  squared  distances  between 
(*o  , ^J# )  and  Sj  ( y4' ,  s* )  : 

CS0i  )*•  (.*•-*£)*■*  (l  »  1,..*,  N)  . 

2.5)  Determination  of  elements  3^j  (i*l, 2 , . . .  ,N;  j»l 

of  matrix  A: 

a^t  s  jlilL.  -  **■'  <. N) 

tsoU>0)a  cs0<:)2 

a  *•-»>•  .  *•_-**»«  U- 1, ,  N). 

(£*o% 

Step  3)  Normal  equations 

3.1)  Determine  matrix  AW  (a  matrix  2  x  N). 

3.2)  Determine  matrix  ATWA  (a  matrix  2x2). 

3.3)  Determine  matrix  (*rW*)  ‘  ’*  (a  matrix  2x2) 

as  indicated  in  step  3.3  of  subsection  II. A. 3. a. 

3.4)  Determine  matrix  /iVl  (a  matrix  2  x  1). 

3.5)  Finally,  determine 

X  -  (  ATW  A  r‘  (  /AT  W'  L  )  . 

Step  4)  First  adjusted  values 


As  indicated  in  Step  4  of  subsection  II. A. 3. a. 


Step  5)  2nd  iteration 

As  indicated  on  Step  5  of  subsection  II. A. 3. a. 
Step  6)  Next  iterations 

As  indicated  on  Step  6  of  subsection  II. A. 3. a. 


C.  FIX  DETERMINATION  BY  TWO  RANGE  DISTANCES  AND  ONE  AZIMUTH 


This  problem  illustrates  how  to  deal  with  observations 
of  different  kinds  (distances  and  angles) .  The  procedures 
for  obtaining  the  residuals  and  the  weight  matrix  are  more 
complex. 

1.  Solution  for  Two  Range  Distances  and  One  Azimuth  from 
3  Different  Stations. 

Given  a  positioning  problem  as  diagramed  in  FIG  II-7, 
in  which: 

R,  -  is  the  observed  range  distance  from  station  #1 

B  l  -  is  the  observed  range  distance  from  station  #2 

/A  -  is  the  observed  azimuth  from  station  #3 
~  are  the  grid  coordinates  of  station  #1 

"  are  the  grid  coordinates  of  station  #2 

C*,H>  -  are  the  grid  coordinates  of  station  #3 


or 

-  is 

the 

standard 

error 

of 

K 

(in  meters) 

<rx 

-  is 

the 

standard 

error 

of 

(in  meters) 

«* 

-  is 

the 

standard 

error 

of 

A 

(in  degrees) 

the  grid  coordinates  of  vessel's  position  P  (x,y)  will  be 
determined. 

Step  1)  Formulation  of  observation  equations 

1.1)  The  analytical  expression  for  the  range  distance 
between  station  i  (i=l,2)  and  vessel's  position  P(xy)  is 
given  by 

rL  *  [  (x- *i.  )*  *  (3  -y  i  )v  ]  ^  s  F  (*,})  (i>  1,2) 


The  function  F(xy)  must  be  expressed  in  a  Taylor's 
series  around  an  "initial  position"  ,  P  ,  whose  coordinates 
are  defined  as  xQ  and  yQ.  Evaluating  the  zero  and  first  order 
terms  of  the  series,  the  following  expression  is  obtained: 

«  [(Yo-  XtjV  4- 


Xo  -  X  c _ 


.  + 


_ 

[CXo-Xt)Z+ 


Then,  designating  by  s^Q  the  distance  from  station  i  (i=l,2) 
to  "initial  point"  PQ  (X0'Y0)  t  the  following  expression  is 
obtained : 

The  observation  equations  are  given  by 

V^-.  i km. a*  +  P*-’3'. .  4?  _  (Pc-  Sio  )  fC=  [,*.) 

Sic  Sic 


where  Rj.  is  the  observed  range  distance. 

In  this  result  it  should  be  noted  that  the  residuals, 
Vj_  (i=l,2),  are  expressed  in  meters. 

1.2)  The  analytical  expression  for  the  azimuth 
between  station  3  and  P(x,y)  is  given  by 

AB  (radians)  =  tan”'  ZLzjQ 

30* 
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Therefore,  the  observation  equation  is  expressed  as 


tan~ 1  *-  _  A  -  Vs 

a  -  TS 

where  A  is  the  observed  azimuth  angle.  In  that  result, it 
should  be  noted  that  is  expressed  in  radians.  Therefore 

it  will  be  necessary  to  obtain  expressed  in  meters. 

(See  FIG  II-5) . 


hjerth 


/  % 


k  A  v 


✓  \ 

;  > 


FIG  II-5 :  CONVERTING  ANGULAR  RESIDUAL  INTO 


METRICAL  RESIDUAL 


From  the  FIG  II-5  is  concluded  that 
(^(meters)  ■  sin  "Vj  (radians)  x  distance  between  P  and 


i-l'4 


V,  ( -  Sir  I"  tin"'  ?***  -  A  +  C3-31) 


Expressing  the  function  Vj  (j$y)  as  a  Taylor's 
series  around  the  "initial  point"  ( *o,yo )  ,  and  taking 
only  the  zero  and  first  order  terms,  the  following  is  obtained 

\JL  _  -  f  f  •»«**  a!  f* .  ..  * 


3z  Sir? 


in  tan'1  -  a]  .  lx0-x ,)  *  (]<,-  it  ^ 

I  *-*  J  L  J 


2*- 

-  * 

Xo  - 

X» 

i 

Xo  - 

X» 

j*- 

*3 

Xp  * 

/3 

33 

Xo- 

Xs 

2)o  "  0^5 


_ „  , 

[(x*-X*iS  (3oOS>*  ]/2> 


- 2^25 _ _  lZ\X  + 

J  (3»^5->Z]  Z  J 

-  ^  * 

J  [" l**-*i?  *  (jo -33)*]^ 


- f*— rniW 

3.- J»  J4, 


Designating  by  S30  and  A*  3  0  the  distance  and  azimuth 
between  station  S3  and  "initial  point”  P0  (vd  ,y0)  , 

then 

>-  r 


S30  s  +  {30  -  }»)*! 


Therefore,  the  observation  equation  may  be  written  as 


s  Sin  $30  +" 

«  oos(A« 3o-  A)-  +  sm  tA^so-A) .  *21*?  .AX  * 

S3o  ^30 

•W(A*j, -/*)•  X»^2  +*» C/l23o  -/»). 

-$3o 

1.3)  Finally, the  observation  equations  in  matrix 
notation  are  expressed  as 

AX  -  L  --  V 

where  the  elements  3,_;j  and  of  matrices  A  and  L 

are  given  by 

cl  II  —  (X©  “  /  S  \q  ^IZ5  /  &\0 

a  (Xo  -  *z)  /  -Sj  a  /  Sio 

a5l*  cos  (A*3o-A)  .  <.  s  11?  (<4*3*-^ .  *Li£! 

S30  S30 

a3z ,  cos  (^130  -  A).  +s;«  [A*30-A). 

$3c  5 3o 

L 1  s  R  |  -  5  ,<>  It  a  Ri  -  52o 

L-.  -  S«*  (  A)  .  S30 


Step  2)  Determination  of  weight  matrix  W 

The  standard  errors  <37  and  Gz  of  range  observations 
are  expressed  in  meters;  the  standard  error  O3  of  the  ob¬ 
served  azimuth  angle  is  expressed  in  degrees.  Therefore,  it 
will  be  necessary  to  obtain  (J3  expressed  in  meters.  (See 
FIG  II-6) . 


FIG  II-6 :  CONVERTING  ANGULAR  STANDARD  DEVIATION 
INTO  METRICAL  STANDARD  DEVIATION 


From  FIG  II-6  is  concluded  that 

O3  (meters)  s  SirrC^  C  degrees)  ^30  . 

Having  obtained  O3  expressed  in  meters,  the  pro¬ 


cedure  for  obtaining  the  weight  matrix  W  is  as  indicated  in 
Step  1  of  subsection  II.A.3.a. 


Step  3)  Normal  equations 


Forming  the  normal  equations,  the  adjusted  values 
for  A  X  and  are  given  by 

X  *  U’WA)'1  UTW  L)  . 

Step  4)  With  the  values  A*  and  A y  a  new  "initial  point" 

Po  (  Is  obtained j 

>  i)o  :  Uj  ♦  4  j  . 

Step  5)  For  the  new  coordinates  (X^jo  )  of  "initial  point", 
the  value  of  ^  (in  meters)  is  recomputed  and  the  weight 
matrix  W  readjusted. 

Step  6)  The  procedure  will  be  repeated,  in  an  iterative  way, 
until  the  increments  AX  and  become  vanishingly  small, 

or,  in  practical  terms,  converging  to  within  a  specified 
tolerance. 

Then,  the  most  probable  values  for  the  coordinates 
(xy)  will  coincide  with  those  obtained  for  the  last  "initial 
point" . 

2.  Numerical  Example 

Referring  to  FIG  II-7,  the  grid  coordinates  (U.T.M.) 
of  shore  stations  are 


COORDINATES 

LUCES  (#1) 

MB 4  (#2) 

MUSSEL  (#3) 

X 

595,794.5 

603,425.2 

597,967.8 

y 

4,055,042.7 

4,053,917.2 

4,053,453.2 

AkRTH 


FIG  II-7 :  FIX  BY  TWO  RANGE  DISTANCES  AND  ONE  AZIMUTH 
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The  following  angle  has  been  measured  by  theodolite: 

P  -  MUSSEL  -  LUCES  *  99®. 360 

For  illustrative  purposes,  the  instrument  is  assigned 
a  standard  error  0"^  =  0*.024. 

The  following  distances  were  measured: 

t 

P  -  LUCES  -  5233  m 
P  -  MB 4  =  3515  m 

Their  standard  errors  are  assumed  to  be  =  Gi  s  1°  m 
(fictitious  values) 

Step  1)  The  azimuth  between  MUSSEL  and  LUCES  is  given  by 

Al3t  *  tanV  -  3*6*  l8l 

3i-  Is 

Therefore , the  following  data  are  available: 


5233  m 

07  a 

l  O  m 

R  2  3 

35  IS  m 

* 

\0  tn 

A  * 

S7* 

C°OZ4 

Step  2)  Formulation  of  observation  equations 

2.1)  Determination  of  first  "initial  point" 

The  first  "initial  point"  will  be  the  point 
determined  by  range  distances  Rj  and  R^,  (for  which  the 
azimuth  from  station  3  is  closer  to  A) .  Therefore,  the  point 
?0  (  *o,3o)  will  satisfy  the  following  system  of  equations: 


(x«  -Xi)  *  “  Ri 

Ua-Xa)*  «  (>t0  if*  *a  . 


Introducing  numerical  values,  the  following 


solution  sets  for  the  above  equations  are  obtained  : 


Jx„,  s  600 }  361.1 

*  V56>320<5 


and 


Jyot3  6oO;2&0-\ 

l^ot  »  kjOS2j3  <i7.  6  # 


The  azimuth  from  station  #3  to  (J^,  1)  and  ^oz,loz^' 

respectively,  are  obtained:  =  45?  24  and  A^3az  * 

115 *.5  .  Therefore,  the  valid  solution  is  the  one  corres¬ 
ponding  to  ,  i.e., 

J  *0=  *o  1  =  600,867.2 
l  3o=  *oi  =  4,056,328.0 

2.2)  Determination  of  azimuth  between  station  3  and 


Po  (*Oo} : 


A*3o  *  tan'1 

3.-33 


Then, 


A* 30  -  4  =  -  d?  297  . 

2.3)  Determining  distances  between  stations  and 

i  K 

5  s  [(  *.-**>a+  C  3,- Jo?  J  -  5Z33-0 

-52o  [(Vi-  V.)%  (3l  -  Jo)1  j14  r  35  15-0 

^  1'/ 

^30  S  £ (^3  "  X«)  +  ( ^3  ~  3o)  J  a  ^  33  .  O  . 


T 


2.4)  Therefore,  the  elements  3 Cj  and  Li,  of  matrices 
A  and  L  will  be 

/  5|  o  s  0.  96  3  3  6  8  Sjja  (ja-CJi}  /  ^  to  =  0-2^5  bIH 

321S  (Xo-XO/Sio  739  3aaa(>-32)/^2o  ^°^8586i 

33|3  <^$(4*30  *>--33  ^  Sot  *kz£  a0.1ool<oo 

53a  Sj© 

3*1:  C£>*(Ai,t  -A)  .  *3-*o  46itxC^3o-A)‘  ~  —  z -0.713*755 

S*o  ^3a 

L|  -  -  5  ,0  2  o 

^2  -  ^2  "  ^2a  =  o 

Ldr  -  5<n  (  Alj0  -  /4  ).  630-  2J* 


2.5)  The  observation  equations  in  matrix  notation 
may  be  written  as 


0.969  3  68  0.  ZWf  614  L  x 

-  0.127  7  38  0.68S  86/ 

0.100  4oa  -0.7<3  75  S  J 
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Step  3)  Normal  equations 


3.1)  Determination  of  weight  matrix  W: 


Then 


SJ7 *  \o  m 
CJs  Id  m 

O3  =*>  ^3  ( Sa*  «  Sin  (d?8  4) 

*  1.710  m 


J  /  <JT  =•  O-  <S7  1 
X  •  0.  O  t 
I  /  <J3'i  r  0.3  4  . 

Setting  the  least  weight  equal  to  one,  it  will  be  obtained 


that 

UJ,  r  » 

UO*  S  | 

OJj  -  34 


or 


o  o 
»  o 
O  34 


3.2)  The  solution  of  normal  equations  is  given 


by 


X=  c/VA)  '  I  a'  WL). 

3.2.1)  Determination  of  A  W 


atw  „ 


0.16^36  8  -0-72773* 
1  o.24S  6/4  0.4B5Q6  l 


23-  80  <5  00 
-24.  24767o 


3.2.2)  Determination 


atwa. 


16. 146 
.-11256 


A  W  A  i 

-  17.25  6 
17.  651  .  . 


3.2.3)  Determination  of 


(ATW*  )*’  : 


0.66  2^5 
0.660  23 


0.  66o  26 
O.  6<»4  2  7 


3.2.4)  Determination  of 


(AWL)  : 


t/WL)3 


5  o4.  ©  15 
-513.  625 

•»  mm 


3.2.5)  Finally, 


Step  4)  First  adjusted  values 

With  the  values  =5.1  and  =  “  23.8  a  new 
"initial  point"  is  obtained  : 

X0  =  600,867.2  +  5.1  s  600,872.3 
>  =  4,056,328.0  -  23.8  =  4,056,304.2 
Step  5)  With  the  new  values  for  the  "initial  point"  the 
procedure  indicated  in  steps  2.2,  2.3,  2.4,  2.5,  3  and  4 
repeated  and  a  "closer"  initial  point  is  obtained. 


is 


Step  6)  That  procedure  must  be  repeated,  in  an  iterative  way, 
until  the  increments  and  become  vanishingly  small, 

or,  in  practical  terms,  converge  to  within  a  specified 
tolerance.  Then,  the  last  "initial  point"  obtained  will 
coincide  with  the  most  probable  position  for  P. 

These  computations  may  be  compared  with  those  shown 
in  the  computer  output  section  on  page  150.  Differences 
in  the  results  are  due  to  the  fact  that  the  calculations 
illustrated  on  the  preceeding  pages  were  only  carried  out  for 
one  iteration. 

3.  Solution  for  the  General  case 

The  solution  will  be  presented  in  such  a  way  that 
easily  can  be  implemented  by  an  algorithm  satisfying  a 
modular  design. 


a.  Two  cases  will  be  considered: 

1st  case)  Two  range  distances  and  one  azimuth  from  three 
stations  (See  FIG.  II-8) 


FI'j  TI-8:  FIX  FROM  3  STATIONS 


Designate  by  S{  and  the  stations  from  which  range 

distances  are  observed;  the  station  from  which  an  azimuth 
is  observed  will  be  designated  by  Sj. 

2nd  case)  Two  range  distances  and  one  azimuth  from  just  two 
stations.  (See  FIG  II-9) 


FIG  II-9:  FIX  FROM  2  STATIONS 

Designate  by  Si  the  station  from  whch  an  azimuth  and  a 
range  distance  were  measured  (  S^concides  with  Sjj )  ;  the 
remaining  station  will  be  designated  by  Sj. 
b.  Given 

(*l,yi)  -  grid  coordinates  of  station  S^ 

(x2,y2)  -  grid  coordinates  of  station  S2 

(X3,yi)  -  grid  coordinates  of  station  S3 

(in  the  2nd  case,  they  coincide  with  the  coordinates 

(X1  /Yl )  of  S]_ ) 

-  range  distance  from  station  S-|_ 

R2  -  range  distance  from  station  S2 

A  -  azimuth  from  station  S-, 


v  i  -  standard  error  of  R]_  (in  meters) 

0 z  -  standard  error  of  R2  (in  meters) 

(J3  -  standard  error  of  A  (in  degrees) 

the  coordinates  of  vessel's  position  P(x,y)  will 
be  determined. 

Step  1)  Formulation  of  observation  equations 

1.1)  Determination  of  first  "initial  point"  {X0^0) 
The  first  "initial  point"  will  be: 

a)  one  of  the  intersection  points  of  circum¬ 
ferences  centered  at  S\  and  with  radius  ranges  of  R, 

and  R2,  respectively. 

b)  the  intersection  point  which  lies  closer  to 
the  azimuth  line  through  S3  . 


Therefore,  the  following  equations  must  be 


satisfied: 


(X,-  to)  t  (  -  R 

(  X2  -  *«)  t  (  3 1  -  Jo)  -  Rt  . 


_  1  _  l  X  i  X  X. 

E  s  R,  -  K  +  32  -  j.  «■  **  -  *  . 


Then, 


E  2  (  3,  -  3i) 


2  tx*-  x.) 


1.1.1)  For  Xi  1*  X\  proceed  as  follows: 


Let 

'  Xx  -  *i 

E2  =•  _  z  x t  r  3'“ 3. z  \  —  L  3  »  t 

(Xi-x',)‘l  ‘  Xi-X.  1 

Ev(— — f-  -  R,% 

v  2(Xz-x,)  •  X2-X, 

EA  =  (  £*)*-  4(e.)  (e,)  . 


Then, 


-  ( E*)  1  7  Eo 

Z  6, 


a)  If  E^<o/  then  the  two  circumfer¬ 
ences  don't  intersect;  therefore,  choose  point  Q  as  first 
"Initial  point" (See  FIG  11-10) . 


^  5<») 

i\ 


S2lXa^x) 


FIG  11-10 :  RANGE  DISTANCES  NOT  INTERSECTING 
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te: 


i  t 


Jo  =  t 


(  y2  -  x  t ) _ 

R  i  (  -  3» ) _ 

[(xl-v,)l+  (3l-3.f  ]/l 


b)  If  E4  *  0  ,  then  the  two  circumferences 
are  tangent.  Therefore 


*.  -  h 


*  E, 


*  O  =  . . .  +  . . . -  jo 

2  ( /z-y.)  **-*. 


c)  If  E4 > o  ,  then  the  two  circumferences 
intersect  at  2  points;  therefore,  the  following  intersection 
points  are  obtained: 


.  [-  ft  We4  ]  /  2  f 


V  E  3.-  2ft  ^ 

2(X*  -*,  ) 


3«, 

]/z  e, 

^  L 

=>'-  31  3at 

2  (Xj  -  X,  ) 

X*  -  *1 

and 


'Soi.  s  ^2  -  >f  ^4  }  /  ^  E  i 


c. 1)  If  X3  =  Xff  |  and  ^3  -  1,  then 

<  X o  -  Xoi 

*•  Jo  -  y  •  *  • 

c. 2)  If  *3  -  XoZ  and  }3=  2*t,  then 


J  X  a  =  Xot 

*  • 

c.3)  otherwise,  determine  azimuths  between 


station  63  (X3,2Ji)  and  iXoL'Xt)  (for 


Criterion: 


D  If  *  ^3  and  >  X*  ,  then 

A23<>c  *  n*/ *• . 

II)  If  J0i  =^3  and  JTdt  <  X3  ,then 

At-zd  *  ?»/t. 

For  ^oi.  ^  ^3  designate  by  o<3t, 

the  solution  given  by  a  computer  of 

^stah'  Xoi-X3  (L-1,1), 
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a)  If  F*o  ,  the  two  circumferences 
do  not  intersect  or  are  tangent; 


then 

*o  -  *i  . 

b)  If  the  two  circumferences 

intersect  at  points  (xaf  and 

?  then 

/x0l  r  X,  4  V  F 

t  -  3* 

and 

(  x©i-  *i  — 

IjoZ*?* 

b .  1 )  If  *3  =  X*  i  and  ^3  =  2©  i  ,  then 

X© 8  Xoi 

b.  2)  If  X3  =  X©2.  and  ^3  *  Joz.  ,  then 

X0 -  Xa 1 

b.3)  Otherwise,  determine  azimuths  A2  301 

and  A 2  joz  between  station 
S3  (x 3,  3>3  )  and  (x»;  ,*•<;  )  (i=l,2) 
using  criterion  presented  in  step  1.1.1 

Having  determined  Ai^0i  and  Ai 1 307.  , 
determine  which  one  is  closer  to  A  . 

I)  If  o\  -  A^ioz  *  A, then 

the  solution  is  undetermined. 


r 


i 

*»» 

R 

* 

** 

,1 

i 


j 

i 

- 

* 

! 


mi"  iw  'V". 


■71  *>}**. 


II)  If  (Aiiifl-A)L  =  (^23oz'AK 

then  choose 

X0  -  *1  . 

III)  If  (A^a-A)'  >  (^23©2 -^)r, 

then  choose 

Xa  -  . 

IV)  If 

then  choose 

X0 =  <oi. 

1.2)  Determining  the  azimuth  A  between 
station  $3  (*3;^3)  and  (X<>;>©) 

Criterion: 

I )  If  "jo  =  ^3  and  X©  >  X3  ,  then 

A*  3d  -  IT' /  *1  . 

II)  If  ?e>=?3  and  X<,<x3)then 

A^3c,  *  3ir  /  £  . 

For  ^  5*  designate  by  0*3  the 
solution  given  by  a  computer  of 

<*3  tan  1  ft-** 

3« 03 
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Ill)  If  and  X0  Xj  jthen 

IV)  If  <©  and  X*<  X3  jthen  A?j0  =  ^3 <2^ 

V)  Otherwise ,  •  ^3*=  ^3  +  ^**  • 

1.3)  Determining  distances  between  stations  and  t 

S  \o  i  C  (*, -**>%  IVfVo? 

S20  ^  C  frt- *•)  t  Oz-3®^1  3  1 

S3*  5  [  Ui-x*f  +  bjo^1  ]  . 

1.4)  Determining  elements  of  matrix  A  , 

a,.  .  *« .-  *1  a,t. 

5  is  s“> 

a*  .  y«„tja.  a„.  2irii 

83 1  j  Cos  (/U  30  *  A).  *5**  (^3o“ 

S50 

a3**  c^sCA^jd-  /0.  3e^ 

S3*  *3« 

1.5)  Determining  elements  of  matrix  L  , 

l|  3  ^  J  "  U  a  &Z  -  $Zo 

1*3  s  Sift  (A  -  A*  3#  )  .  5  3^ 

Step  2)  Solution  of  normal  equations 

2.1)  Obtain  weight  matrix  W. 

2.1.1)  Determine  standard  error  of  observed  azimuth 
angle  expressed  in  meters, 

O3  (meters)  s  sin  C3  (radians)  x  5 30 


i - :  n 

1  AD-A123  917  CALCULATION  OF  HVDROGRAPHIC  POSITION  DATA  BV  LEAST  ' 

P  SQUARES  ADJUSTMENT CU)  NAVAL  POSTGRADUATE  SCHOOL 

NONTEREV  CA  F  CASTRO  E  SILVA  JUN  82 

I  UNCLASSIFIED  F/G  12/1 

flBHl 

u 

1 

■ 

— — 

_ 

2.1.2)  With  U3  expressed  in  meters ,  the  procedure 
for  obtaining  the  weight  matrix  W  is  as 
indicated  in  Step  1  of  subsection  ll.A.3.a 

2.1.3)  Finally/  (Tj  is  again  expressed  in  radians; 

(J~\  (  r  aAians)  s  Sin  £  /  SgoJ 

2.2)  Determine  matrix  ArVl. 

2.3)  Determine  matrix  fiTWA  . 

2.4)  Determine  matrix  (At  WA)~  I 

2.5)  Determine  matrix  L  . 

2.6)  Finally/  determine  (AT*WL)# 

Step  3)  First  adjusted  values 

With  the  values  AX  and  A-^)  obtain  new  "initial 
point"  Pq  ^o); 

U-  Xe+  A* 
l  y p  -  do+  a ^  • 

Step  4)  2nd  iteration 

For  obtaining  a  "closer"  initial  point, repeat  steps 
1.2,  1.3,  1.4,  1.5,  2  and  3. 

Step  5)  Next  iterations 

Repeat  step  4  until  AX  and  become  vanishingly 

small,  or,  in  practical  terms,  converge  to  within  a  specified 
tolerance.  Then,  the  adjusted  values  for  x  and  y  will 
coincide  with  the  coordinates  of  the  last  "initial  point" 


obtained. 


Ill  .  RESULTS  AND  CONCLUSIONS 

A.  RESULTS 

From  the  general  case  solutions  developed  for  the  selected 
positioning  methods,  algorithms  were  written  in  a  structured 
programming  format.  All  algorithms  are  presented  in  appendix 
I. 

These  modular  algorithms  were  translated  into  Fortran 
language  for  implementation  on  the  NPS  computer,  an  IBM  3033. 
Program  listings  are  provided  in  the  Computer  Programs  Section 
beginning  on  page  151. 

Data  sets  given  in  each  Numerical  Example  section  were 
input  into  the  corresponding  computer  program,  and  the  output 
of  each  run  is  given  in  the  Computer  Output  Section  starting 
on  page  148. 

Additionally,  the  programs  were  tested  using  several  fic¬ 
titious  data  sets  to  insure  their  performance  in  handling  the 
various  initial  conditions  which  were  modeled  for  each  fix¬ 
ing  method. 

In  applying  these  programs  to  real  positioning  data  the 
following  points  should  be  considered: 

1.  The  presence  of  blunders  and  systematic  errors  in 
the  observations  will  be  reflected  in  the  dimensions  of  the 
error  ellipse.  If  all  blunders  are  removed  by  careful 
editing  and  all  systematic  errors  are  eliminated  by  modeling 
or  calibration,  then  the  size  of  the  error  ellipse  will 


represent  the  positioning  error  due  to  net  geometry  and 
random  errors. 

2.  When  the  information  about  the  standard  errors  of  the 
observations  is  reliable  (for  example,  determined  by  field 
calibration  procedures) ,  then  the  estimates  obtained  for  the 
standard  deviations  of  the  observed  values  will  be  close  to 
the  a  priori  values  (see  example  of  fix  by  3  azimuth  angles 
in  Computer  Output  section  on  page  148). 

3.  When  no  a  priori  values  are  given  for  the  standard 
deviations  of  the  observations  (it  is  assumed  that  the 
observations  are  equally  weighted) ,  then  the  application  of 
the  least  squares  method  will  provide  estimates  of  instrument 
(or  observation)  accuracy  (see  example  in  computer  output 
section,  on  page  149). 

4.  When  the  correlation  coefficient  is  close  to  one, 
the  error  ellipse  becomes  flatter  approaching  a  straight 
line  (see  example  of  fix  by  two  range  distances  and  one 
azimuth  angle  in  the  Computer  Output  section  on  page  150) . 

5.  When  the  correlation  coefficient  is  negative,  the 
major  axis  of  the  error  ellipse  runs  through  the  2nd  and 
4th  quadrants.  Thus,  the  angle  from  the  x-axis  to  the  major 
axis  measured  counterclockwise  lies  between  90°  and  180°.  If 
the  correlation  coefficient  is  positive,  the  major  axis  runs 
through  the  1st  and  3rd  quadrants,  and  the  angle  from  x-axis 
to  the  major  axis  measured  counterclockwise  is  between  0°  and 


B.  CONCLUSIONS 


The  most  significant  result  of  this  thesis  is  that  well 
documented  programs  are  now  available  which  can  be  used  for 
the  analysis  of  hydrographic  positioning  data.  These  programs 
may  be  employed  to  process  and  analyze  hydrographic  survey 
data  that  have  been  collected  using  one  of  the  three  position¬ 
ing  methods  discussed.  Ideally,  such  software  should  be 
adapted  to  run  in  a  mini  computer  aboard  a  survey  vessel  or 
launch.  This  capability  would  allow  "real  time"  analysis  of 
positioning  accuracy. 

In  addition  to  processing  actual  survey  data,  the  programs 
may  assist  in  survey  planning.  By  scaling  observations  from 
existing  charts  of  a  survey  area,  sample  data  sets  may  be 
formed  to  test  net  geometry.  This  information  can  be  used 
to  establish  the  best  location  for  shore  control  stations. 

The  programs  are  written  in  modular  form  so  that  they 
may  be  adapted  for  use  by  other  types  of  positioning  systems. 
The  significant  differences  between  all  the  programs  lie  in 
the  modules  dealing  with  the  computation  of  the  "initial 
point"  and  formulation  of  the  observation  equations. 

It  should  be  noted  that  the  accuracy  of  the  geodetic 
control  stations  has  not  been  specifically  considered  in 
these  formulations.  However,  any  survey  error  in  the  station 
coordinates  will  be  reflected  in  the  dimensions  of  the  error 
ellipse  of  the  adjusted  hydrographic  position. 


All  of  the  programs  were  developed  using  a  plane 
coordinate  system  model.  Thus,  they  are  primarily  applicable 
to  nearshore  hydrographic  positioning  problems.  Application 
to  offshore  hydrography  would  require  a  geodetic  coordinate 
system  model  based  on  a  selected  spheroidal  datum  surface. 
Obviously , the  use  of  a  geodetic  coordinate  system  would  yield 
more  complex  analytical  expressions  relating  the  unknowns. 
But,  once  these  were  obtained  and  linearized,  then  the  pro¬ 
cedures  for  computing  adjusted  survey  coordinates  and  the 
statistical  values  defining  their  precision  are  identical  to 
those  developed  in  this  thesis. 

Whether  the  existing  programs  are  used  in  their  current 
form  or  modified  to  accomodate  other  variables,  one  final 
point  should  be  made.  The  most  significant  contribution  of 
the  least  squares  method  to  hydrographic  position  adjustment 
is  its  ability  to  quantify  errors  statistically.  When 
programs  are  operated  aboard  the  survey  vessel  in  "real  time 
relative  accuracy  achieved  with  conventional  survey  methods 
is  elevated  to  absolute  accuracy  if  redundant  observations 
are  made  and  adjusted  using  least  squares. 

Monitoring  the  size  and  orientation  of  the  error  ellipse 
alerts  the  user  to  the  presence  of  gross  blunders  and  inord¬ 
inately  large  systematic  errors.  The  need  for  electronic 
positioning  system  calibration  can  be  realistically  evaluated 
and  calibration  may  be  performed  on  an  as  needed  basis.  With 


sufficient  redundant  observations,  electronic  positioning 
systems  can,  in  fact,  become  self  calibrating. 

As  the  trends  in  electronic  and  computer  technology 
continue  to  decrease  the  cost  of  collecting  and  processing 
redundant  observations,  conventional  two  LOP's  survey 
positioning  will  be  relegated  to  the  historical  equivalent 
of  lead  line  hydrography. 
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APPENDIX  A.  LEAST  SQUARES  PRINCIPLE  AND  NORMAL  DISTRIBUTION 
When  measuring  a  parameter,  the  outcomes  of  that  experi¬ 
ment  can  be  considered  as  values  assumed  by  a  random  variable 
following  a  normal  distribution.  For  a  random  variable  X 
following  a  normal  distribution,  the  value  most  likely  to 
occur  is  its  mean  ft*  .  The  true  value,  from  a  deterministic 
point  of  view,  of  an  observed  parameter  is  ,  in  a 
stocastical  sense  the  mean  of  the  random  variable  associated 
with  the  experiment.  Therefore,  when  using  the  least  squares 
technique  for  the  adjustment  of  a  redundant  number  of  obser¬ 
vations,  not  only  a  set  of  "consistent"  values  are  ob¬ 
tained  but  also  the  most  probable  values  for  the  means  of 
the  random  variables  considered.  Therefore,  the  adjusted 
values  are  also  the  best  estimates  for  the  "true"  values  of 
the  parameters  considered. 

1.  Normal  distribution 

The  density  function  associated  with  a  random  variable 
X  following  a  normal  distribution  is  expressed  by  (see 
FIG  A-l) . 

-  [(X-/.,)1  /l  < T?  ] 

i*)  *  — 1 -  e 

Ox  VTtr 
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FIG  A- 3:  RESIDUALS 


If  the  vessel  is  situated  at  0,  then  the  means  of  the  random 
variables  X,  Y  and  Z,  associated  with  the  range  distances 
AO,  BO  and  CO,  are  at  distances  V^,V 2  and  V3  from,  respec¬ 
tively,  observed  values  x,  y  and  z.  (See  FIGS  A-3  and  A-4) . 


LO 


FIG  A-4:  RESIDUALS  AND  NORMAL  DISTRIBUTION 


The  probability  of  occurrence  of  a  set  of  observed  values 
■{x  ^  x}  is  given  by 


j?l  y  x}  )  s  R  Px  Uii})  .  Pa  l  {*))  . 

If  the  observations  are  equally  weighted  then  the  standard 
deviations  <£  ,  fly  and  fl^  have  the  same  value,  say,  . 
Then,  from  (A-l)  ,  it  will  be  obtained 


Pi\x*  x))*(\cl<r>/ZZ)  e 


-  [(*-/•*)*  {?->»*>*♦ 


C^-2) 


Recalling  that 

V,  =  y-r.  -  * 
Vi*  M  ~  ' 1 
\4  =/**-»  , 


it  will  be  obtained  from  (A-2) 


P  (  |  X  (j  )  r  CoMSTAAiT  X 


0*3) 


The  means  A  ,  and  ^L»|  are  values  such  that  the  observed 

values  *,  y  and  *  have  a  probability  as  high  as  can  be 
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APPENDIX  B.  LEAST  SQUARES  PRINCIPLE  FOR  WEIGHTED  OBSERVATIONS 
1.  If  an  observed  value  Xt‘  has  the  weight  ,  then  the 

observed  value  "X(  is  worth  as  much  as  observed  values-^ 
with  weight  equal  to  unity. 

A  usual  criterion  for  establish  weights  is  to  consider 
the  weights  inversely  proportional  to  the  squared  standard 
deviations ,  i .  e .  , 


UJ, 


Then,  given  a  set  of  observed  values  XL*  (i=0,l,2, .  . . )  ,  un¬ 
equally  precise,  with  standard  deviations  ,  if  the  least 
precise  value  x0  is  considered  to  have  a  weight  equal  to 
unity  ^0=1)  ,  the  different  weights  UL),;  will  satisfy 

(f—  Z  •  X 

-  V—  l  1*0,1, 

<rr‘ 


2.  A  set  of  observed  values  x^ ,  x2 • • • • >  xn ,  with  respective 
weights  equal  to  a),  ,  , . .  .  ,  n  >  is  equal  to  a  set  of 

"l  values  equal  to  ,  uJj  values  equal  to  , .  . .  ,  and 
values  equal  to  ,  all  with  unity  weight.  Therefore,  the 

sum  of  the  squared  residuals  will  be  given  by 


2  Z  x"2- 

t  /x,  -  *.  )  +  "t  K /Xx  -  •  •  •  4  ( ./■*„  -*•«), 

and  the  basic  least  squares  principle  will  be  expressed  as 

b  n 


ZLj  i  M/.  *  2  V*  - 

L;i  /  u  i.  *  I 


mm  imi)  n 
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APPENDIX  C.  NORMAL  EQUATION  IN  ALGEBRAIC  NOTATION 
Below  are  the  observation  equations  in  algebraic  notation: 


Vi  *  3 1|  *  i 

+  atl 

4  3  l  m  *  -n  “ 

L,  =  o 

,  Vzsaai*, 

t  3az 

*x  4  .  .  . 

4  3lm  zm  - 

L  l  s  O 

.  ^n-3Blx, 

4  3 n2 

xz  .  •  . 

+  3  n  m  x  m  ” 

Ln  =  ° 

The  unknown  values  X, ,  x x  that  satisfy  the  basic 

least  squares  principle 


minimun 


L  »  » 


are  those  that  satisfy  the  following  expression: 

LMi  (3  ii  xi  4  3ri  xi  +  ••  •  +  3im  Lj)  + 

W-j  at, ^xa.  4  ...  4  dxtm  3c pn~  )  4 

_ _ _ _ _ _ .4 

^nUr,,  *«43«Z  ^  4  ...  4  3„mxm  -Ln)  *»•••*!»»)* 

=  TtiLnCrnurj  . 

The  values  minimizing  F(x^,x2»***  are  those 

such  that 


Considering  that 


^JL  a  2  [  3  c  l  }  y,  . .  •  +  2  [  3  Cl  3  Xm  "  2.  [u>t-  3t*t  1 t 1  =  ® 

3X, 

IE  *  1  [U1taw»3ii]yi4*-*  44  [u>i.ata  a<lm  ]  3li  ^u]=  0 

ax,  ” 


£E.  =  Z  [“‘i.ai1aim ]*,+  ...  -tzK  a;^  J*.-1  [“>;  3i-*  Li  1  -  0 

d*m 

the  following  normal  equations  are  obtained  : 


[“>i  au‘  . 

*i  t.  .  aim!  Km  m[wC  aU  u] 

[mi  at,  aii  . 

X,  «■  .  .  .  +[u>L3it  acm]xm  Li! 

j wc  au,  *i  J*,  + « •  -  *  [<*>c  atn?  !*•»  ’  LuiSim  ^J=  ° 


APPENDIX  D.  NORMAL  EQUATIONS  IN  MATRIX  NOTATION 
Taking  the  observation  equations  in  its  matrix  notation, 


AX-  L»  V, 

the  unknown  vector  X  satisfying  the  basic  least  squares 
principle  , 

VT  W  V  -  miLnonun  , 

is  the  one  such  that 

(.ax-  litw  (/» y-  l)  = 

XT/1TW  A  X  -  1  XT/ATW  L  *  lT  W  L.  =  minimun  . 

The  vector  X  satisfying  the  above  expression  will  be  such 
that 

^ _  (xVwAX-2XVWLt  LtWU=o  .  (D-1 

c>X 

Considering  that 

]L  (XtAtWAX)  =  2  XTArWA 

ax 

(  xtatw  l)  =  lt  wj 
ax 

JL  (  lt  W  l  )  =  o  , 

0X 

it  will  be  obtained  from  (D-l)  that 

1  xVwA  -  2  l_TWr  A  =  o  . 


1 


Therefore,  the  normal  equations  are  of  the  form 


(AtW  A)  X  *  ATW  L, 

and  the  solution  will  be  given  by 

X  •  iaVa)'1  (AtV  L  ). 
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APPENDIX  E.  A  COMPUTATIONAL  CHECK  FOR  THE  LEAST  SQUARES 
Adjustment  Technique 
By  taking  the  normal  equations 

[Ar w  A)  X  =■  AT W  L  (  E  -  1) 

and  recalling  that 

A  *  =  V+  L, 

then  the  following  can  be  obtained : 

(AtWa)X-/wV.  AtWL.  (£-2) 

Therefore,  from  (E-l)  and  (E-2)  it  is  obtained  that 

ATWV'=0.  (£-3) 

Equation  (E-3)  provides  a  check  on  the  computations  for  least 
squares  adjustment. 
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APPENDIX  F.  THE  CONTROVERSIAL  CRITERION  FOR  ASSIGNING  WEIGHTS 
1.  The  usual  criterion  for  assigning  weights  is  stated  as: 

a)  the  weights  are  inversely  proportional  to  the  squared 
standard  deviations,  i.e., 

tu;  f  (F-i; 

where  K  is  an  arbitrary  constant? 

b)  the  least  precise  observation  has  the  unity  weight, 

1. e. , 

<Jj0  So  =  K  =>  K  =  C 1  )  .  or 

K  -  So  (F-l} 

where  Sd  is  the  standard  deviation  of  least  precise  obser¬ 
vation. 

2.  For  the  moment,  consider  only  equation  (F-l) .  The 
influence  of  the  value  assigned  to  k  on  the  computations 
for  obtaining  the  adjusted  values  and  standard  deviations  of 
adjusted  values  will  be  determined.  From  eq  (F-l)  it  will  be 
obtained  that 

-  KT  /  Su  . 

Then,  the  weight  matrix  W  will  be 
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and  the  trace  TR  {W  )3  1C  f  $*  t  ..  .  ♦  K  /  5«  a  K  Tft  (W#)  . 


a)  Computing  adjusted  values, 

X  --  (  atv/o-'  (atwl>  » 

*  1'/k)(atw'a)‘!  kiatw'u  *  iaVa)"1  (atw*l). 

Therefore,  it  is  concluded  that,  for  the  adjusted 
values  X,  the  value  assigned  to  K  is,  in  fact,  arbitrary. 

b)  Computing  standard  deviations  of  adjusted  values. 


4.  -  4  A 


where 


is  an  element  of  matrix 


Recalling  that 

Q.  (AtWaj'=  ( i/k j  MTWAr‘a  (i/k)  a1, 

then 

=  C  •  /  K)  q'it  . 

Next,  by  considering 

5es  /vTwv  .  Jk  .  I  vrW  , 

V  TRCNA/j-m  V  K  TR(w')-  m 

it  follows  that 

SXi-  A'ii 


VT  W  1  V 
K  TRW')- m 


Therefore,  as  should  be  expected,  the  standard  deviations  of 
adjusted  values  are  affected  by  the  value  assigned  to  K. 


3.  In  fact,  eq  (F-2)  imposes  a  constraint  on  the  K  value: 


To  illustrate  the  consequences  of  accepting  that  kind  of 
constraint,  suppose  that,  given  100  observations,  99  are 
equally  precise  and  one  is  less  precise,  say,  with  a  standard 
deviation  400  times  greater  than  the  standard  deviation  of  the 
remaining  99  observations.  Then,  according  to  the  usual 
criterion,  the  least  precise  observation  has  the  unity  weight 
and  the  other  99  observations  have  the  weight  20.  That  dis¬ 
tribution  of  weights  does  not  seem  "good,"  and  it  is  the 
author's  opinion  that  there  should  be  a  better  constraint 
minimizing  the  disturbances  introduced  by  the  assignment  of 
different  weights. 
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Consider  L\  ,  *-2.  and  L3  as  values  assumed,  respectively, 
by  independent  random  variables  L/  ,  L2  andLj*  Since 
the  observations  were  equally  weighted,  then  L,,La  and 
L3  present  the  same  standard  deviation,  say,  SQ.  Then, 

VAR(X)  =  H  oic  XAR  (  Li)  -  L*4  J  5a 


R  (  y)  =  £  A*  VAR  (  U)  =  [p>1  1  50t  ((9-2 


Co'JA  R 


a 

(X,?)  -  21  <*ty8;  V/*R(  Li)  =■  [=1/9  ]  5.T  . 


.1.-1 

2.  Determining  the  matrix  Q  =  (A  W  /I  )  ,  the  result  is 


MM -[3b]1 

[3']Lbl]-[abf 

-[3b] 

Ca*l 

-  [a  b] 1  (3’][b?]  -  (abj 


Since 


U]. 

03  .  -  „  „ 

&‘][b']-[ab]» 

L*1)- 

[a\)0>‘]-[abr  " 

1 

V 

O' 

1 

D 

k 

C,X][bI]-[ab]  * 

it  may  be  concluded  from  eq.  (G-2)  that 


|  s  S0  .  »• 

^  Sy  ;  So  .  ^ 

Sx;y  ?  5e  •  . 


4.  Finally,  the  correlation  coefficient  between  random 
variables  x  and  y  is  given  by 


P:  _ 

’  5 x  «  5* 
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Then,  the  probability  of  occurrence  of  the  vessel's  position  in 
a  small  area  (44*.  44^)  aroundf(x,y)  will  be  given  by 


(see  FIG  H-l) 

P  C  x=*tA*,  y=^±^j)=  /  /  d3ax 

“jC-a* 


2.  Now,  it  will  be  determined  in  what  kind  of  line  points  with 
equal  probability  of  occurrence  are  situated. 

These  points  will  present  the  same  value  f  for  the 
density  function.  Therefore,  letting 


1C » 


x. 


K2  =  [-  in  ( 2  tr  s[K\  §  )  (2  Ki )  J  /  (Tx 

and  iC3=  -  K*  <TX  2  ^  then 

(TylX*-  4.  (T^  >cy+-lQ'jl+-k'3-0  (H-1) 

That  quadratic  equation  in  x  and  y  represents  a  conic;  the 

existence  of  the  xy-term  indicates  that  the  conic  is  rotated 

out  of  its  standard  position. 

a.  Before  determining  what  kind  of  conic  equation  (H-l) 

represents,  check  if  points  {TV)o)  and  (o;<F^)  are  both 

over  the  same  contour  line  for  a  constant  density  function  . 

Inserting  point  {%t}o)  into  (H-l)  results  in 

K3=  -<T^<Ty  \ 

Insertina  point  (O.Jv  )  into  (H-l)  results  in 

h  =  -r;  <T*\ 

Therefore,  the  points  (<3y;0)  and  ( oJ<Jy  )  are  over  the 
same  contour  line  (corresponding  to  ^3=  —  ffx  )  • 


b.  The  analytical  expression  for  the  specific  contour 


line  containing  points 

(<J*  j  0  )  and 

(o,r^  ) 

is 

i 

Oy  x  - 1  <S^  *  3 

+  cr;^rt-(rx'L 

_  2 

{p;  =  o 

(H- 

Considering 

> 

it 

3 

a 

Cr  <rS 

0 

II 

LU 

3  = 

0  =  o 

p= 

it  is  concluded  that  the  discriminant  is  less  than  or  equal 
to  zero,  i.e. , 

B*-  4  A  c  ^  o  . 

7 

Therefore,  if  d-*i  Ac<o  then  equation  (H-2)  represents  an 

t 

ellipse;  if  3 -Mc  =  o  then  it  will  represent  a  straight 
line  (a  degenerate  ellipse  corresponding  to  a  perfect 
correlation  between  random  variables  x  and  y) . 

3.  The  equation  of  the  error  ellipse  in  standard  position: 
Consider 


(  AtW  A  Q 


^  i  ^3 
^3  <U 


Recalling  that 


i- 

t- 

t*' 

4 


4 


then  equation  (H-2)  is  equivalent  to 

S  *  -  z  *3  4  «5ol  -  *  .  (  H-3) 

Consider 

A-u  Cx%  E=  « 

B  X  -1<|»  0.0  F*  -q,  lx  S.* 

If  the  x— axis  is  rotated  an  angle  'tfo  such  that 

cot  Igi  r  diC  a  q.-qi  (H-4) 

B  <*3 

then  the  equation  of  the  ellipse  (H-3)  in  its  standard 
position  is 


(  So’  \ 


t 


q.<u  s.1 


(H-5) 


M3  $*«»&*  +  S»  C«al*  3* 


•) 


After  some  algebraic  and  trigonometric  manipulation,  the 
following  expression  is  obtained  from  (H-5) 


-  D 


(  SoZ  ) 

*U  ♦  0 


1 


(H-6) 


S 
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where 


D- [  t«|i  -  q*  )*+  4  ]  .  (H-7) 

4.  If  the  positive  value  of  D  satisfying  eq.  (H-7)  is 
choosen,  then  the  semi-major  axis  of  the  error  ellipse  is 
positioned  along  the  "new"  x-axis.  Therefore,  from  the 

A 

two  solutions  Of,  r  of  and  ot,r  of*  9o  satisfying  eq.  (H-4)  the 
valid  one  must  be  choosen  (considering  o(  as  the  smallest 
positive  angle  satisfying  (H-4) )  . 


FIG  H-2 :  ERROR  ELLIPSE 


The  only  way  to  solve  that  ambiguity  is  to  test  either  with 
ol,  or  o^. 


For  that  purpose,  it  is  recommended  to 

a)  obtain  the  point  of  intersection  (either  P, 


or  t*  ) 


of  the  line  y  *  xi an*  with  the  ellipse  before  rotation 
(expressed  by  eq  (H-3)); 

b)  determine  the  distance  _d  between  the  origin  and 
(or  ); 

c)  if  d  =  semi-major  axis  a  ,  then  the  major  axis  makes 

an  angle  (measured  counterclockwise)  with  the  "old" 

x-axis;  if  not,  then  the  angle  will  be  0/2  =0/1+90  , 

that  is,  the  major  axis  runs  through  2nd  and  4th  quadrant. 
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APPENDIX  I  -  ALGORITHMS 


MODULE  I 

ALGORITHM  FIX-BY-N-AZIMUTHS 
INPUT  N 

PI-3.141592653589793 
OUTPUT 'NUMBER  OF  STATIONS3', N 
DO  FOR  I— 1  TO  N 

INPUT  TABLE-INPUT (1,1), TABLE-INPUT (1,2), TABLE-INPUT (1,4) 
OUTPUT ' ST# ' , I , ' EAST= ' , TABLE-INPUT (1,1),' NORT3 ' , 
TABLE-INPUT (I, 2), 'ST  ERROR3' , TABLE-INPUT (I, 4) 

END  DO 

DO  FOR  1-1  TO  N 

INPUT  TABLE-INPUT (I, 3) 

END  DO 

DO  WHILE  TABLE-INPUT (1,3) *400.0 
OUTPUT 'OBSERVED  AZIMUTHS' 

DO  FOR  I«*-l  TO  N 

OUTPUT 'AZIMUTH  FROM  STATION# ’ , I , ' 3 ' , TABLE-INPUT (I , 3 ) , 

’ DEGREES ' 

ALGORITHM  CONVERSION-DEGREES-RADIANS (TABLE-INPUT (I , 3 ) ) 
TABLE-INPUT (1,3) ^TABLE-INPUT ( I , 3 ) * ( PI / 1 80 . 0 ) 

END  CONVERSION-DEGREES-RADIANS (TABLE-INPUT (1,3)) 

END  DO 

MODULES  2, 3, 4, 5, 6, 7 
DO  FOR  I  — 1  TO  N 

INPUT  TABLE-INPUT (I, 3) 

END  DO 
END  DO 

END  FIX-BY-N-AZIMUTHS 
MODULE  2 

ALGORITHM  WEIGHT -MATRIX  (N , TABLE-INPUT ( I , 4) ) 

ALGORITHM  ZERO  (TABLE -WEIGHT) 

DO  FOR  I— 1  TO  10 
DO  FOR  J— -1  TO  10 

TABLE-WEIGHT ( I , J ) —  0 . 000 
END  DO 
END  DO 

END  ZERO  (TABLE-WEIGHT) 

ALGORITHM  SQUARE  (N , TABLE-INPUT( I , 4) , TABLE-WEIGHT) 

DO  FOR  I- 1  TO  N 

TABLE-WEIGHT (1,1)*- TABLE-INPUT ( I , 4)**2 
END  DO 

END  SQUARE  (TABLE-WEIGHT) 

ALGORITHM  NORMALIZE  (TABLE-WEIGHT) 

GREATEST  -  TABLE-WEIGHT (1,1) 

DO  FOR  I  *-2  TO  N 

IF  TABLE-WEIGHT (1,1)  >  GREATEST  THEN 
GREATEST  «- TABLE-WEIGHT (1,1) 

END  IF 


DO  FOR  I«-l  TO  N 

TABLE-WEIGHT (1,1) ^-GREATEST /TABLE-WEIGHT (1,1) 

END  DO 

END  NORMALIZE  (TABLE -WEIGHT) 

END  WEIGHT -MATRIX  (TABLE-WEIGHT) 

MODULE  3 

ALGORITHM  FIRST-INITIAL-POINT  (TABLE-INPUT ,N) 

ALGORITHM  SELECT-AZIMUTHS (TABLE-INPUT (1,3) ,N) 

I«-2 

DO  WHILE  TANGENT (TABLE-INPUT (1,3) ) =TANG ENT (TABLE - 
INPUT (1,3)) 

Iv-I+1 

IF  2>N  THEN 

OUTPUT’ POSITION  IS  UNDETERMINED  FOR  THAT  DATA  SET' 
PICK  UP  ANOTHER  DATA  SET 
END  IF 
END  DO 

END  SELECT-AZ IMUTH ( I ) 

ALGORITHM  INITIAL-COORDINATES  (TABLE-INPUT , I ) 

IF  TABLE-INPUT ( 1 , 3 ) *0 . 0  OR  TABLE-INPUT (1,3) “PI  THEN 
MK  «- TANGENT ( ( 5 . / 2 . )*PI -TABLE-INPUT (1,3)) 

XO  ♦•  TABLE-INPUT (1,1) 

YO  <*•  TABLE-INPUT (1,2)  +MK*  (XO-TABLE-INPUT  (1,1)) 

ELSE  IF  TABLE-INPUT (I, 3) *0.0  OR  TABLE-INPUT (I , 3 ) =PI  THEN 
MI  ♦•TANGENT ( (5 . /2 . )*PI-TABLE-INPUT( 1 ,3) ) 

XO  4-TABLE-INPUT (1,1) 

YO  ♦•TABLE-INPUT (1,2) +MI* (XO-TABLE-INPUT (1,1)) 

ELSE 

MI  ♦•  TANGENT (  (5  .  /2  .  )*PI-TABLE-INPUT(  1 , 3)  ) 

MK  ♦-  TANGENT ( ( 5 . / 2 . )*PI-TABLE-INPUT (1,3)) 

XO  ♦-  (TABLE-INPUT (1,2) -TABLE-INPUT (1,2) +MI*TABLE- 
INPUT (1,1) -MK*TABLE-INPUT (1,1))/ (MI-MK) 

YO  ♦-TABLE-INPUT (1,2) +MI* (XO-TABLE-INPUT (1,1)) 

END  IF 

END  INITIAL-COORDINATES  (XO,YO) 

END  FIRST-INITIAL-POINT (XO,YO) 

MODULE  4 

ALGORITHM  ITERATIONS  (TOLERANCE) 

DO  UNTIL  TOLERANCE  <1.0 

MODULES  8,9,10,11,12,13,14 
END  DO 

END  ITERATIONS (XO , YO ) 

MODULE  5 

ALGORITHM  FINAL-ADJUSTED-POSITION  (XO,YO) 

OUTPUT’ ADJUSTED  COORDINATES  X* ’ ,XO, ' Y* ' ,YO 
END  FINAL-ADJUSTED-POSITION  (X,Y) 
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MODULE  6 

ALGORITHM  PRECISION  (TABLE-A, TABLE -WEIGHT, TABLE -Q, LI ST-L, 

DELTAX , DELTAY , N ) 

MODULES  15,16,17,18,19 
END  PRECISION (SU , SX, SY , SXY ,RO) 

MODULE  7 

ALGORITHM  ERROR-ELIPSE (TABLE -Q , SU ) 

OUTPUT' ERROR  ELIPSE  SEMI -AXIS  AND  ORIENTATION' 

MODULES  20,  21,  22,  23,  24,  25,  26 
END  ERRO-ELIPSE(SA, SB ,GAMAO) 

MODULE  8 

ALGORITHM  INITIAL-AZIMUTHS (XO , YO , TABLE-INPUT, N) 

DO  FOR  I«-l  TO  N 

IF  YO^T ABLE- INPUT (1,2)  AND  XO  >  TABLE-INPUT (1,1)  THEN 
LIST-AO(I) *“PI/2 . 

ELSE  IF  Y0=TABLE-INPUT(I,2)  AND  XO < TABLE-INPUT (I , 1 )  THEN 
LIST-AO(I)  «-  (3 .0*PI)/2  . 

ELSE  IF  XO=TABLE-INPUT (1,1)  AND  YO  > TABLE-INPUT (I , 2)  THEN 
LIST-AO(I)  4-0.0 

ELSE  IF  XO=TABLE-INPUT (1,1)  AND  YO  <TABLE-INPUT(I , 2 )  THEN 
LIST-AO(I)*-  PI 
ELSE 

ALFA(I)  4- ARC  TANGENT ( (XO-TABLE-INPUT(1 , 1 ) ) / (YO- 
TABLE-INPUT (1,2))) 

IF  ALFA(I)  >0.0  AND  XO  > TABLE-INPUT (I , 1 )  THEN 
LIST-AO(I)*-  ALFA(I) 

ELSE  IF  ALFA(I)  >0.0  AND  XO  <  TABLE-INPUT (1,1)  THEN 
LI  ST -AO  ( I )  ♦-  ALFA  ( I )  +PI 

ELSE  IF  ALFA(I )  <0.0  AND  XO  >  TABLE-INPUT  (1 , 1)  THEN 
LIST-AO(I)  4-  ALFA(I)+PI 
ELSE 

LI ST-AO ( I ) +- ALFA ( I ) + 2 . 0*P I 
END  IF 
END  IF 
END  DO 

END  INITIAL-AZIMUTHS (LIST-AO) 

MODULE  9 

ALGORITHM  MATRIX-L (TABLE-INPUT (1,3), LIST-AO , N ) 

ALGORITHM  ZERO(LIST-L) 

DO  FOR  I«-l  TO  10 
LIST-L(I)  4-0.000 
END  DO 

END  ZERO(LIST-L) 

DO  FOR  I  •- 1  TO  N 

LIST-L( I ) 4- TABLE-INPUT ( I , 3 )-LIST-AO(I ) 

END  DO 

END  MATRIX-L (LIST-L) 
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MODULE  10 

ALGORITHM  INITIAL-SQUARED-DISTANCES (XO,YO,N, TABLE-INPUT) 

DO  FOR  I*-l  TO  N 

LIST-SO(I)  *-(XO-TABLE-INPUT(1 , 1) )**2+(YO-TABLE- 
INPUT(I,2))**2 

END  DO 

END  INITIAL-SQUARED-DISTANCES (LIST-SO) 

MODULE  11 

ALGORITHM  MATRIX-A  (N , TABLE-INPUT ,XO ,YO , LIST-SO) 

ALGORITHM  ZERO  (TABLE-A) 

DO  FOR  14-1  TO  10 
DO  FOR  J  4-1  TO  2 

TABLE-A<  I,  J)  *-0.000 
END  DO 
END  DO 

END  ZERO (TABLE-A) 

ALGORITHM  ELEMENTS (TABLE-A ,N , TABLE-INPUT ,XO ,YO,LIST-SO) 
DO  FOR  I  *-l  TO  N 

TABLE-A (1 , 1 )  4-  ( YO-TABLE-INPUT (1,2)) /LIST-SO ( I ) 

END  DO 

DO  FOR  I*-l  to  N 

TABLE-A (1,2)  4- (TABLE-INPUT ( 1 , 1 ) -XO ) /LIST-SO ( I ) 

END  DO 

END  ELEMENTS (TABLE-A) 

END  MATRIX-A (TABLE-A) 

MODULE  12 

ALGORITHM  NORMAL-EQUATIONS  (TABLE -A, TABLE-WEIGHT, LI ST-L) 
ALGORITHM  ATW(TABLE-A , TABLE-WEIGHT) 

DO  FOR  I*-l  TO  2 

DO  FOR  J*-l  TO  10 

TABLE-ATW  ( I ,  J  )  4-  0 . 00 
DO  FOR  K-  1  TO  10 

TABLE-ATW  ( I ,  J  )  *-  TABLE-ATW  (I  ,J)  +TABLE  - 
A(K,I)*TABLE-WEIGHT(K,J) 

END  DO 
END  DO 
END  DO 

END  ATW( TABLE-ATW) 

ALGORITHM  ATWA  (TABLE-ATW, TABLE-A) 

DO  FOR  I  ♦-  1  TO  2 
DO  FOR  J4-1  TO  2 

TABLE-ATWA  ( I ,  J )  *-0.00 
DO  FOR  K*-l  TO  10 

TABLE-ATWA (I , J)  TABLE-ATWA (I , J ) +TABLE- 
ATW(I,K)*TABLE-A(K,J) 

END  DO 
END  DO 
END  DO 

END  ATWA  (TABLE-ATWA) 


ALGORITHM  INVERT -ATWA (TABLE -ATWA) 

BETA «— TABLE-ATWA ( 1 , 2 )**2-TABLE-ATWA( 1 , 1 ) *T AB LE -ATWA ( 2 , 2 ) 
TABLE-Q( 1 , 1 )*-  -TABLE -ATWA ( 2 , 2 ) /BETA 
TABLE-Q (1,2)*- TABLE -ATWA ( 1,2) /BETA 
TABLE-Q (2,1)*-  TABLE-Q (1,2) 

TABLE-Q (2 , 2 ) *- -TABLE-ATWA ( 1 , 1 ) /BETA 
END  IN VERT -ATWA (TABLE-Q) 

ALGORITHM  ATWL(TABLE-ATW,LIST-L) 

DO  FOR  I*  1  TO  2 

LIST-ATWL(I) *-0.0 
DO  FOR  K*-l  TO  10 

LIST-ATWL(I) * LIST-ATWL(I)+TABLE-ATW(I ,K)*LIST-L(K) 
END  DO 
END  DO 

END  ATWL (LI ST-ATWL ) 

ALGORITHM  AD JUSTED- INCREMENTS ( TABLE-Q , LI ST-ATWL ) 

DELTAX  «-  TABLE-Q ( 1 , 1 )*LI ST-ATWL ( 1 ) +TABLE- 
Q( 1 , 2)*LIST-ATWL(2 ) 

DELTAY  *- TABLE-Q ( 2 , 1 )*LIST-ATWL( 1 ) +TABLE- 
Q ( 2 , 2 )*LIST-ATWL ( 2 ) 

END  ADJUSTED-INCREMENTS (DELTAX, DELTA  Y) 

END  NORMAL-EQUATIONS (DELTAX, DELTA  Y) 

MODULE  13 

ALGORITHM  NEW-INITIAL-POINT  (XO ,YO .DELTAX, DELTAY) 

XO-XO+DELTAX 
YO  *-YO+DELTAY 

END  NEW-INITIAL-POINT(XO.YO) 

MODULE  14 

ALGORITHM  TOLERANCE (DELTAX, DELTAY) 

TOLERANCE*-  DELTAX**2+DELTAY**2 

END  TOLERANCE (TOLERANCE) 

MODULE  15 

ALGORITHM  RESIDUALS (TABLE-A , LIST-L , DELTAX , DELTAY , N ) 

LIST-X(l)*- DELTAX 
LIST-X(  2)*- DELTAY 

ALGORITHM  LIST-AX (TABLE-A , LIST-X ,N) 

DO  FOR  1*1  TO  N 
LIST-AX(I)*-0.0 
DO  FOR  J  *  1  TO  2 

LIST-AX(I) *  LIST-AX(I)+TABLE-A(I ,J)*LIST-X(J) 

END  DO 
END  DO 

END  LIST-AX (LIST-AX) 


ALGORITHM  LIST-V(LIST-AX,LIST-L ,N) 

DO  FOR  I  ^ 1  TO  N 

LIST-V (I)  *- LIST-AX ( I ) -LIST-L ( I ) 

END  DO 

END  LIST-V(LIST-V) 

END  RESIDUALS (LIST-V) 

MODULE  16 

ALGORITHM  ST-DEVIATION-OF-UNIT-WEIGHT-OBS ( LIST-V, TABLE-WEIGHT, N) 
ALGORITHM  LIST-VTW(LIST-V , TABLE-WEIGHT , N ) 

DO  FOR  I  *- 1  TO  N 

LIST-VTW (I) ♦- LIST-V ( I ) *TABLE-WEIGHT (1,1) 

END  DO 

END  LIST-VTW (LIST-VTW) 

ALGORITHM  VTWV (LIST-VTW, LIST-V) 

VTWV  <r  0 . 0 
DO  FOR  I*-l  TO  N 

VTWV  *- VTWV+LIST-VTW ( I )*LIST-V ( I ) 

END  DO 

END  VTWV (VTWV) 

ALGORITHM  TRACE (TABLE-WEIGHT) 

TRACE  **0.0 

DO  FOR  I  *-l  TO  N 

TRACE  *-TRACE+TABLE-WEIGHT(I , I) 

END  DO 

END  TRACE (TRACE) 

ALGORITHM  SU ( VTWV , TRACE) 

CHARLIE  <r (VTWV/ (TRACE-2 . 0 ) ) 

SU  *■  SQRT  ( CHARLIE ) 

END  SU(SU) 

END  ST-DEVIATION-OF-UNIT-WEIGHT-OBS ( SU ) 

MODULE  17 

ALGORITHM  ST-DEVIATION-OF-EACH-OBS ( SU , TABLE-WEIGHT ) 

OUTPUT' PRECISION  OF  OBSERVATIONS' 

DO  FOR  I<-1  TO  N 

S  <-(SU/ SQRT (TABLE-WEIGHT (I , I))) *(180. 0/PI) 

OUTPUT 'ST  DEVIATION  OF  OBS ', I ,'=', S DEGREES ' 

END  DO 

END  ST-DEVIATION-OF-EACH-OBS 
MODULE  18 

ALGORITHM  ST-DEVIATIONS-AND-COVARIANCE-OF  X-AND-Y ( SU , TABLE-Q ) 

SX  *-  SU*SQRT (TABLE-Q (1,1)) 

SY  *-  SU*SQRT  (TABLE-Q  (2,2)) 

SXY  *-  (SU**2)*TABLE-Q( 1,2) 

OUTPUT 'SX=' ,SX, ’SY=' ,SY, 'SXY=' ,SXY 
END  ST-DEVIATIONS-AND-COVARIANCE-OF-X-AND-Y ( SX , SY , SXY ) 


MODULE  20 

ALGORITHM  D(TABLE-Q) 

D  4-SQRT ( (TABLE-Q (1,1) -TABLE-Q ( 2 , 2 ) )**2+4 . 0* (TABLE- 
Q(l,2)**2)) 

END  D(D) 

MODULE  21 

ALGORITHM  SEMI-MAJOR-AXIS ( SU , TABLE-Q) 

SA4-  SU*SQRT  ( 2  .  Q*TABLE-Q  ( 1 , 1 )*TABLE-Q  (  2 , 2  )/<  TABLE-Q  (1,1)  + 
TABLE-Q (2, 2 )-D)) 

OUTPUT 'SEMI -MAJOR  AXIS  SA=',SA 
END  SEMI-MAJOR-AXIS (SA) 

MODULE  22 

ALGORITHM  SEMI-MINOR-AXIS ( SU , TABLE-Q ) 

SB  «-SU*SQRT  (  2 . 0*TABLE-Q  (1,1  )*TABLE-Q  (2,2  )/( TABLE-Q  (1,1)  + 
TABLE-Q ( 2 , 2 ) +D ) ) 

OUTPUT ' SEMI -MINOR  AXIS  SB®* ,SB 
END  SEMI-MINOR-AXIS (SB) 

MODULE  23 

ALGORITHM  GAMA (TABLE-Q) 

IF  TABLE-Q (1,1) =TABLE-Q (2,2)  THEN 
GAMA«-PI/4.0 
ELSE 

OMEGA  4- ARC  TANGENT ( 2 . 0*TABLE-Q (1,2)/ (TABLE-Q (1,1)- 
TABLE-Q (2,2))) 

IF  OMEGA  >0.0  THEN 
GAMA «— OMEGA/  2 . 0 
ELSE 

GAMA  <r  (OMEGA+PI )  /  2 . 0 
END  IF 
END  IF 

END  GAMA (GAMA) 

MODULE  24 

ALGORITHM  INTERSECTION ( SU , TABLE-Q , GAMA ) 

X10  «-( SU**2 )*TABLE-Q ( 1 , 1)*TABLE-Q (2,2) 

XI 1*—TABLE-Q ( 2 , 2 ) -2 . 0*TABLE-Q ( 1 , 2)*TANGENT (GAMA) + 
(TANGENT (GAMA*^2)*TABLE-Q( 1 ,1) 

X1^~X10/X11 

Y1 «- XI* (TANGENT (GAMA)** 2) 

END  INTERSECTION  (Xl.Yl) 

MODULE  25 

ALGORITHM  AVERAGE (S A, SB) 

AVER ♦- ( (SA+SB) /2 . Q)**2 
END  AVERAGE (AVER) 
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MODULE  26 

ALGORITHM  SELECTION (AVER , XI ,Y1) 

Dl «“  Xl+Yl 
IF  Dl >  AVER  THEN 
GAMAO  ♦*  GAMA 
ELSE 

GAMAO  «-GAMA+PI/2 .0 
END  IF 

GAMAO  GAMAO*  ( 180 . 0  /  PI ) 

OUTPUT'ANGLE  FROM  X-AXIS  TO  SA  ANTICLOCKWISE® ' , GAMAO 
END  SELECTION (GAMAO) 

MODULE  19 

ALGORITHM  CORRELATION-COEFFICIENT ( SX , SY , SXY) 
R04-SXY/(SX*SY) 

OUTPUT ' CORRELATION  COEFFICIENT  RO=',RO 
END  CORRELATION-COEFFICIENT (RO) 
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MODULE  30 

ALGORITHM  FIX-BY-N-SEXTANT-ANGLES 
PI  ♦*3.141592653589793 
INPUT  N 

OUTPUT 'NUMBER  OF  SEXTANT  ANGLES =',N 
M  **N+1 

DO  FOR  I  4-1  TO  N 

INPUT  TABLE-INPUT (1,1), TABLE-INPUT (1,2), TABLE- 
INPUT(I ,  4) 

OUTPUT 'ST#' ,1, 'EAST*' , TABLE-INPUT (I, IX' NORTH*’ , TABLE- 
INPUT  ( I , 2 ) , ' ST  ERROR* ’ , TABLE-INPUT (1,4) 

END  DO 

INPUT  TABLE-INPUT (M , 1 ) , TABLE-INPUT (M , 2 ) 

OUTPUT' ST#' ,M, 'EAST*' , TABLE-INPUT (M, 1) , 'NORTH*' , TABLE- 
INPUT  (M,  2) 

DO  FOR  14-1  TO  N 

INPUT  TABLE-INPUT (I, 3) 

END  DO 

DO  WHILE  TABLE-INPUT (1,3) *400.0 
OUTPUT 'OBSERVED  SEXTANT  ANGLES’ 

DO  FOR  14-1  TO  N 
J4-I+1 

OUTPUT 'SEXTANT  ANGLE  BETWEEN  ST#', I,' AND  ST#’,J, 

' * ’ , TABLE-INPUT (1,3),' DEGREES ’ 

ALGORITHM  CONVERSION-DEGREES-RADIANS (TABLE-INPUT (1,3)) 
TABLE-INPUT (1,3)  TABLE-INPUT ( I , 3 )* ( PI / 1 80 . 0 ) 

END  CONVERSION-DEGREES-RADIANS (TABLE-INPUT (I , 3 ) ) 

END  DO 

MODULES  2,31,32,5,6,7 
DO  FOR  14-1  TO  N 

INPUT  TABLE-INPUT(I , 3 ) 

END  DO 
END  DO 

END  FIX-BY-N-SEXTANT-ANGLES 
MODULE  31 

ALGORITHM  FIRST-INITIAL-POINT-FOR-FIX-BY-N-SEXTANT- 
ANGLES (TABLE-INPUT ) 

MODULES  33,34,35,36 

END  FIRST-INITIAL-POINT-FOR-FIX-BY-N-SEX'’’ANT-ANGLES  (XO  .  YO  ) 
MODULE  32 

ALGORITHM  ITERATIONS (TOLERANCE ) 

DO  UNTIL  TOLERANCE  <1.000 

MODULES  37,38,39,40,12,13,14 
END  DO 

END  ITERATIONS (XO.YO) 


MODULE  33 

ALGORITHM  SELECT-SEXTANT -ANGLES (TABLE-INPUT ) 

J*-l 

DO  UNTIL  FRAC  1  f  FRAC  2 
J«-J  +  l 

IF  J>M  THEN 

OUTPUT' SOLUTION  UNDETERMINED  FOR  THAT  DATA  SET' 

PICK  UP  ANOTHER  DATA  SET 
END  IF 
I  — J-l 
K— J  +  l 

ANGUL  TABLE-INPUT (1,3) +TABLE-INPUT ( J , 3 ) 

FRAC  1  *■  COSINE (ANGUL)*SQRT ( ( (TABLE-INPUT (1 ,  1 )- 
TABLE-INPUT (J , 1 ) )** 2+ (TABLE-INPUT (I , 2 ) -TABLE- 
INPUT  (J , 2) )**2)* ( (TABLE-INPUT (K , 1 ) -TABLE-INPUT (J , 1 ) )**2+ 
(TABLE-INPUT (K , 2 ) -TABLE-INPUT (J , 2 ) )**2 ) ) 

FRAC  2  *-( TABLE-INPUT (I , I ) -TABLE-INPUT (J , 1 ) )* (TABLE- 
INPUT  (J , 1 ) -TABLE- INPUT (K , 1 ) ) + (TABLE-INPUT (1,2) 
-TABLE-INPUT (J ,2) )* (TABLE-INPUT ( J , 2 ) -TABLE-INPUT (K , 2 ) ) 
END  DO 

END  SELECT-SEXTANT-ANGLES ( I, J,K) 

MODULE  34 

ALGORITHM  INTERCHANGE-DATA (TABLE-INPUT , I , J ,K) 

STORE ( 1 ) -  TABLE-INPUT (1,1) 

STORE ( 2 )—  TABLE-INPUT (1,2) 

STORE ( 3 ) —  TABLE-INPUT (1,3) 

STORE (4)4- TABLE-INPUT (2,1) 

STORE (5)4-  TABLE-INPUT (2,2) 

STORE ( 6 )  .-TABLE-INPUT (2,3) 

STORE ( 7 ) .-TABLE-INPUT (3,1) 

STORE (8)*— TABLE-INPUT (3,2) 

STORE ( 9 )  .—TABLE-INPUT (1,1) 

STORE ( 10  )—  TABLE-INPUT (1,2) 

STORE ( 1 1 )—  TABLE-INPUT (1,3) 

STORE( 12 )—■  TABLE-INPUT (J , 1 ) 

STORE (13)4- TABLE-INPUT (J , 2 ) 

STORE (14).-  TABLE-INPUT (J ,3) 

STORE  ( 15  ) —TABLE-INPUT  (K ,  1 ) 

STORE (16).- TABLE-INPUT (K , 2 ) 

TABLE-INPUT (1,1)  -STORE(9) 

TABLE-INPUT (1,2)4-  STORE ( 10 ) 

TABLE-INPUT (1,3 )  .-STORE (11) 

TABLE-INPUT  (2,1)  —  STORE  (12) 

TABLE-INPUT  (2,2)  —  STORE  (13) 

TABLE-INPUT (2, 3)— STORE( 14) 

TABLE-INPUT (3,1)—  STORE (15) 

TABLE-INPUT  (3,2)  —  STORE  (16) 

END  INTERCHANGE-DATA (TABLE-INPUT , STORE ) 
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ALGORITHM  INITIAL-COORDINATES (TABLE-INPUT ) 

IF  TABLE-INPUT  (1 ,3)?‘(PI/2 .0)  AND  TABLE-INPUT(2 ,3)^(PI/ (2 .0) 

THEN 

AB  ^-TANGENT (TABLE-INPUT (1,3)) 

BA  «- TANGENT (TABLE- INPUT (2,3)) 

E  — (TABLE-INPUT (2 , 1)-TABLE-INPUT(1 , 1) )/AB+ (TABLE- 
INPUT  (2,1) -TABLE-INPUT (3,1))/ BA+TABLE-INPUT (3,2)- 
TABLE, INPUT (1,2) 

F  —  (TABLE-INPUT (2,2) )-TABLE-INPUT (1,2)) /AB+ (TABLE-INPUT (2.2)- 
TABLE-INPUT (3,2))/ BA+TABLE-INPUT (1,1)- 
TABLE-INPUT (3,1) 

G  — (TABLE-INPUT (1,2 )*TABLE-INPUT (2,1) -TABLE-INPUT (2,2)* 

TABLE-INPUT (1,1)) /AB+( TABLE-INPUT (2,1 )*TABLE-INPUT (3,2)- 
TABLE-INPUT(3 , 1 )*TABLE-INPUT(2 , 2) ) /BA+TABLE- 
INPUT  (2,1 )*TABLE-INPUT (3,1) +TABLE-INPUT (2,2 )*TABLE- 
INPUT (3,2) -TABLE-INPUT (1,1 )*TABLE-INPUT (2,1)- 
TABLE-INPUT ( 1 , 2 )*TABLE-INPUT( 2 , 2 ) 

IF  F=0 .0  THEN 
DAO  —  G/E 
Y01  —  DAO 
Y02  -  DAO 
U— AB 

R  —TABLE-INPUT (1,2) -TABLE-INPUT (2,2) -AB* (TABLE-INPUT (1,1)  + 
TABLE-INPUT (2,1)) 

SAL-  AB* (DAO**2-DAO* (TABLE-INPUT (1,2) +TABLE-INPUT (2,2))+ 
TABLE-INPUT (1,1 )*TABLE-INPUT (2,1) +TABLE-INPUT (1,2)* 
TABLE-INPUT (2,2)) +DAO* (TABLE-INPUT (2,1) -TABLE-INPUT (1,1))+ 
TABLE-INPUT (1,1 )*TABLE-INPUT (2,2) -TABLE-INPUT (2,1)* 
TABLE-INPUT (1,2) 

DISC  — SQRT(R**2-4.0*U*SAL) 

X01  — (-R+DISC)/(2.0*U) 

X02  —  (-R-DISC)/(2.0*U) 

ELSE  IF  E=0 . 0  THEN 
H  — (-G/F) 

X01  — H 
X02  — H 
U  AB 

R  —  TABLE-INPUT ( 2,1) -TABLE-INPUT ( 1 , 1 )-AB* (TABLE-INPUT ( 1,2)+ 
TABLE-INPUT (2, 2)) 

SAL- AB* (H**2-H* (TABLE-INPUT ( 1 , 1 ) +TABLE-INPUT ( 2 , 1 ))  + 

TABLE-INPUT (1,1 )*TABLE-INPUT (2,1) +TABLE-INPUT (1,2)* 
TABLE-INPUT (2,2)) +H* ( TABLE-INPUT (1,2) -TABLE-INPUT (2,2))+ 
TABLE-INPUT ( 1 , 1 )*TABLE-INPUT(2 , 2) -TABLE-INPUT (2 , 1 )* 
TABLE-INPUT (1,2) 

DISC  — SQRT(R**2-4 . 0*U*SAL) 

Y01  —  ( -R+DISC ) / ( 2 . 0*U ) 

Y02  (-R-DISC)/(2.0*U) 

ELSE 

C—F/E 


136 


DAO  *“G/E 

U  4-  AB*(C**2+1.0) 

R  *-  AB* ( 2 . 0*C*DAO-C* (TABLE-INPUT (1,2) +TABLE-INPUT (2,2)) 
-TABLE-INPUT (1,1) -TABLE-INPUT (2,1)) +C* (TABLE- 
INPUT  (2,1) -TABLE-INPUT (1,1)) +TABLE-INPUT (1,2)- 
TABLE-INPUT (2,2) 

SAL  — AB* (DAO**2-DAO* (TABLE-INPUT (1,2) + TABLE-INPUT (2,2))  + 
TABLE-INPUT (1,1 )*TABLE-INPUT (2,1) +TABLE-INPUT (1,2)* 
TABLE-INPUT (2,2)) +DAO* (TABLE-INPUT (2,1) -TABLE- 
INPUT  (1,1)) +TABLE-INPUT (1,1 ) *TABLE- INPUT (2,2)- 
TABLE-INPUT (2,1) *TABLE-INPUT (1,2) 

DISC  — SQRT(R**2-4 . 0*U*SAL) 

XOl — ( -R+DISC ) / ( 2 . 0*U ) 

X02  — (-R-DISC)/(2.0*U) 

YOl— C*X01+DAO 
Y02— C*X02+DAO 
END  IF 

ELSE  IF  TABLE-INPUT ( 1 ,3)=(PI/2.Q)  AND  TABLE-INPUT (2, 3)* 

(PI/2.0)  THEN 

BA  +-  TANGENT (TABLE-INPUT (2,3)) 

E  —  BA* (TABLE-INPUT (3, 2 ) -TABLE-INPUT ( 1,2))+ 

TABLE-INPUT (2,1) -TABLE-INPUT (3,1) 

F  BA* (TABLE-INPUT (1,1 )-TABLE-INPUT(3 , 1 ) )+TABLE-lNPUT (2,2)- 

TABLE-INPUT (3,2) 

G  —  BA* (TABLE-INPUT (2,2 )*T ABLE-INPUT (3,2) +TABLE-INPUT (2,1)* 
TABLE-INPUT (3,1) -TABLE-INPUT (1,2) *TABLE-INPUT (2,2)- 
T ABLE-INPUT (1,1 )*TABLE-INPUT (2,1))+ 

TABLE-INPUT (2,1 )*TABLE-INPUT (3,2) -TABLE-INPUT (3,1)* 
TABLE-INPUT (2, 2) 

IF  F=0 . 0  THEN 
DAO+-G/E 
YOl-DAO 
Y02-DA0 

R  —  -TABLE-INPUT (1,1) -TABLE-INPUT (2,1) 

SAL  DAO** 2 -DAO* (TABLE-INPUT (1,2) +TABLE-INPUT (2,2))  + 
TABLE-INPUT (1,2) *TABLE-INPUT (2,2) +TABLE- 
INPUT (1,1) *TABLE-INPUT (2,1) 

DISC  *-SQRT(R**2-4.Q*SAL) 

X01+-(-R+DISC)/2.0 
X02—  (-R-DISC)/2.0 
ELSE  IF  E=0 .0  THEN 
H  *-(-G/F) 

X01-H 
X02  -H 

R  +- -TABLE-INPUT (1,2) -TABLE-INPUT (2,2) 

SAL  «- H**2-H* (TABLE-INPUT (1.1) +TABLE-INPUT (2.1)  + 
TABLE-INPUT (1,1) )*TABLE-INPUT (2,1) +TABLE- 
INPUT (1,2 )*TABLE- INPUT (2,2) 

DISC«*SQRT(R**2-4 . 0*SAL) 

YOl-  ( -R+DISC) /2 .0 
Y02  *■  (-R-DISC)  / 2 . 0 
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ELSE 

C  — F/E 
DAO  — G/E 
U+-C**2+1.0 

R  —2 . 0*C*DAO-C* (TABLE-INPUT ( 1 , 2 ) +TABLE-INPUT ( 2 ,  2  ) ) 
-TABLE-INPUT (1,1) -TABLE-INPUT (2,1) 

SAL  *■  DAO** 2 -DAO* (TABLE- INPUT (1,2) + TABLE- INPUT ( 2 , 2 ) )  + 
TABLE-INPUT (1,2) *TABLE-INPUT (2,2) +TABLE-INPUT (1,1)* 
TABLE-INPUT (2,1) 

DISC  *“SORT (R**2-4 . 0*U*SAL) 

XOl— (-R+DISC)/(2.0*U) 

X02 — (-R-DISC) / (2 .0*U) 

YO 1  —  C*XQ 1  +DAO 
Y02— C*X02+DAO 
END  IF 

ELSE  IF  TABLE-INPUT (1,3)?* (PI/ 2.0) AND  TABLE-INPUT ( 2 , 3 ) » 
(PI/2.0)  THEN 

AB  TANGENT (TABLE-INPUT (1,3)) 

E  —  AB* (TABLE-INPUT (1,2) -TABLE-INPUT (3,2)) +TABLE- 
INPUT(1,1)-TABLE-INPUT(3,1) 

F  ■*-  AB*  (TABLE-INPUT  (3,1)  -TABLE-INPUT  (1,1))  +TABLE- 
INPUT (1,2) -TABLE-INPUT (2,2) 

G  —AB* (TABLE-INPUT (1,1 )*TABLE-INPUT (2,1) +TABLE-INPUT (1,2)* 
TABLE-INPUT (2,2) -TABLE-INPUT (2,1 )*TABLE-INPUT (3,1)- 
T ABLE-INPUT (2,2 )*TABLE-INPUT (3,2)) +TABLE-INPUT (1,1)* 
TABLE-INPUT (2,2) -TABLE-INPUT (2,1 )*T ABLE-INPUT (1,2) 

IF  F=0.0  THEN 
DAO-G/E 
Y01  —DAO 
Y02  —  DAO 

R  * — TABLE-INPUT (2,1) -TABLE-INPUT (3,1) 

SAL  — DA0**2-DA0* ( TABLE-INPUT (2,2) +TABLE-INPUT (3,2))  + 
TABLE-INPUT (2,2 )*TABLE-INPUT (3,2) +TABLE-INPUT (2,1)* 
TABLE-INPUT (3,1) 

DISC— SQRT(R**2-4.0*SAL) 

XOl  — (-R+DISC)/2.0 
X02—  (-R-DISC)/2.0 
ELSE  IF  E*0 . 0  THEN 
H*~(-G/F) 

XOl  —  H 
X02  *”  H 

R-  -TABLE-INPUT (2,2) -TABLE-INPUT (3,2) 

SAL-  H**2-H*(TABLE-INPUT(2 , l)+TABLE-INPUT(3, 1 ) )+ 

TABLE-INPUT (2,1 )*TABLE-INPUT (3,1) +TABLE-INPUT (2,2)* 
TABLE-INPUT (3, 2) 

DISC  -  SQRT(R**2-4 . 0*SAL) 

Y01 — (-R+DISC) /2 . 0 
Y02  — (-R-DISC) /2 .0 
ELSE 

C  *“F/E 
DAO— G/E 
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U-C**2+1.0 

R  *-  2 . 0*C*DAO-C* (TABLE-INPUT (2,2) +TABLE-INPUT (3,2))- 
T ABLE-INPUT (2,1) -TABLE-INPUT (3,1) 

SAL  «- DAO** 2 -DAO* (TABLE-INPUT (2,2) +TABLE-INPUT (3,2)) +TABLE 
INPUT (2,1 )*TABLE-INPUT (3,1) +TABLE-INPUT (2,2)*  TABLE- 
INPUT  (3, 2) 

DISC*-SQRT(R**2-4.0*U*SAL) 

X01  ♦-  ( -R+DISC ) / ( 2 . 0*U ) 

X02— (-R-DISC)/(2.0*U) 

YOl  +-(C*X01+DAO 
Y02  +-C*X02+DAO 
END  IF 
ELSE 

E  TABLE-INPUT (1,2) -TABLE- INPUT (3,2) 

F  ^-TABLE-INPUT (3,1) -TABLE-INPUT (1,1) 

G  *- TABLE-INPUT (1,1) *TABLE-INPUT (2,1) +TABLE- 

INPUT (1,2 )*TABLE-INPUT (2,2) -TABLE-INPUT (2,2  +TABLE- 
INPUT  (3,2) -TABLE-INPUT (2,1 +TABLE-INPUT (3,1) 

IF  F=0 .0  THEN 
DAO  +-G/E 
YOl  ♦•DAO 
Y02**  DAO 

X01  ♦-TABLE-INPUT (1,1) 

X02*“  TABLE-INPUT (1,1) 

ELSE  IF  E“0.0  THEN 
H  ♦”  (-G/F) 

X01*-H 

X02*-H 

YOl  ♦-TABLE-INPUT (1,2) 

Y02  ♦-TABLE-INPUT (1,2) 

ELSE 

C  *”  F/E 
DAO  ♦G/E 
U-  C**2  +  l  .0 

R  —  2 . 0*C*DAO-C* ( TABLE-INPUT (2,2) +TABLE-INPUT (3,2))- 
TABLE-INPUT ( 2 , 1 )-TABLE-INPUT (3,1) 

SAL  -  DA0**2-DA0*  (TABLE-INPUT  (2,2)  •'■TABLE- 
INPUT  (3,2))  +T  ABLE-INPUT  (2,2)  ^ABLE- 
INPUT  (3,2) +TABLE-INPUT (2,1 +TABLE-INPUT (3,1) 

DISC  *-SQRT(R**2-4.0*U*SAL) 

X01  ♦-(-R+DISC)/ (2 .0*U) 

X02*-(-R-DISC)/ (2 .0*U) 

YOl  *-C*X01+DAO 
Y02*-C*XQ2+DA0 
END  IF 
END  IF 

ALGORITHM  SELECTION(TABLE-INPUT ,X01 ,X02 , YOl ,Y02) 

IF  TABLE-INPUT (1,3 )^(PI/ 2.0)  THEN 

VALOR  1  *-  ( (TABLE-INPUT ( 2 , 1 ) -XO 1 )* (TABLE-INPUT (1,2)- 
Y01 )-(TABLE-INPUT( 1 , 1 )-X01 )* (TABLE-INPUT (2,2)- 


Y01) ) / ( (TABLE-INPUT ( 2 , 2 ) -YOl )* (TABLE-INPUT (1,2)- 
YO 1 ) + (TABLE-INPUT ( 2 , 1 ) -X0 1 )* (TABLE-INPUT (1,1)- 

xoi)) 

VALOR  2  *-( (TABLE-INPUT (2 ,1)-X02)* (TABLE-INPUT (1 ,2)- 
Y02 ) - (TABLE-INPUT ( 1 , 1 ) -X02 ) * ( TABLE-INPUT (2,2)- 
Y02) )/( (TABLE-INPUT (2,2) -YQ2 )* (TABLE-INPUT (1,2)- 
Y02 ) + (TABLE-INPUT (2,1) -X02 )* (TABLE-INPUT (1,1) 
-X02) ) 

MODUL  1*-ABS (TANGENT (TABLE-INPUT (1, 3 )) -VALOR  1) 
MODUL  2  4r-ABS (TANGENT (TABLE-INPUT(  1,3) )-VALOR  2) 

IF  MODUL  1<  MODUL  2  THEN 
X0*-X01 
YO  Y01 
ELSE 

XO  <-X02 
YO  «-Y02 

END  IF 

ELSE  IF  TABLE-INPUT (2 ,3)^(PI/2.0)  THEN 

VALOR  i  4- ( (TABLE-INPUT ( 3 , 1 ) -XO 1 )* (TABLE-INPUT (2,2)“ 

YO 1 ) - (TABLE-INPUT ( 2 , 1 ) -XO 1 )* (TABLE-INPUT (3,2)- 
YOl ) ) / ( (TABLE-INPUT ( 3 , 2 ) -YOl )* (TABLE-INPUT (2,2)- 
YO 1 ) + (TABLE-INPUT ( 3 , 1 ) -XO 1 )*TABLE-INPUT ( 2 , 1 ) “XO 1 ) ) 
VALOR  2  4-  ( (TABLE-INPUT (3 , 1)-X02)* (TABLE-INPUT (2,2)- 
Y02 )- (TABLE-INPUT ( 2 , 1 ) -X02 )* (TABLE-INPUT (3,2)- 
Y02 ) ) / ( (TABLE-INPUT ( 3 , 2 ) -Y02 )* (TABLE-INPUT (2,2)- 
Y02 )  +  ( TABLE-INPUT ( 3 , 1 ) “XO  2 ) * ( TABLE-INPUT ( 2 , 1 ) -XO  2 ) ) 
MODUL  1 *- ABS (TANGENT (TABLE-INPUT (2,3)) -VALOR  1) 

MODUL  2  ABS (TANGENT (TABLE-INPUT (2, 3 ))-VALOR  2) 

IF  MODUL  1 <  MODUL  2  THEN 
X04-X01 
YO  4-YOl 
ELSE 

XO  «-X02 
YO  4-Y02 
END  IF 
ELSE 

VALOR  1  *-(TABLE-INPUT(2,2)-Y01)*(TABLE-INPUT(l,2)- 
YO 1 ) + (TABLE-INPUT ( 2 , 1 ) -XO 1 ) * (TABLE-INPUT ( 1 , 1 ) -XO 1 ) 
VALOR  2  4-  (TABLE-INPUT ( 2 , 2 ) -Y02 )* (TABLE-INPUT (1,2)- 
Y02 ) + (TABLE-INPUT ( 2 , 1 ) -X02 ) * (TABLE-INPUT ( 1 , 1 ) -X02 ) 
MODUL  l«r- ABS  (VALOR  1) 

MODUL  2*-ABS(VALOR  2) 

IF  MODUL  1  <  MODUL  2  THEN 
XO  *-XQl 
YO*- YOl 
ELSE 

XO  4-X02 
YO  *-Y02 
END  IF 
END  IF 
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END  SELECTION  (XO,YO) 

END  INITIAL-COORDINATES (XO,YO) 

MODULE  36 

ALGORITHM  RESTORE-INITIAL-DATA (TABLE-INPUT , STORE) 

TABLE-INPUT  ( 1 , 1 )  *-  STORE  ( 1 ) 

TABLE-INPUT (1,2)4- STORE ( 2 ) 

TABLE-INPUT (1,3)*-  STORE ( 3 ) 

TABLE-INPUT  (2,1)*-  STORE  (  4  ) 

TABLE-INPUT  (2,2)4-  STORE  ( 5 ) 

TABLE-INPUT (2,3)*- STORE ( 6 ) 

TABLE-INPUT (3,1)*-  STORE ( 7 ) 

TABLE-INPUT  (3,2)*-  STORE  (  8  ) 

TABLE-INPUT  (1,1)*-  STORE  ( 9  ) 

TABLE-INPUT (1,2)*-  STORE (10) 

TABLE-INPUT  (1,3)*-  STORE  (11) 

TABLE-INPUT ( J , 1 )  4-  STORE (12) 

TABLE-INPUT(J , 2 ) *-STORE( 13 ) 

TABLE-INPUT (J ,3) — STORE (14) 

TABLE-INPUT (K , 1 ) *- STORE (15) 

TABLE-INPUT (K, 2) *- STORE (16) 

END  RESTORE-INITIAL-DATA (TABLE-INPUT) 

MODUDULE  37 

ALGORITHM  INITIAL-AZIMUTHS (XO , YO , TABLE-INPUT , M) 

DO  FOR  I«-l  TO  M 

IF  Y 0 *TABLE - INPUT (1,2)  AND  XO  > TABLE-INPUT (1,1)  THEN 
LIST-AZ(I)  **(3 .0*PI)/2 .0 

ELSE  IF  YO=TABLE-INPUT(I,2)  AND  XO < TABLE-INPUT (I , 1 )  THEN 
LIST-AZ (I)*- PI/2.0 
ELSE 

LIST-ALFA ( I )  *-  ARC  TANGENT ( (TABLE-INPUT ( 1 , 1 ) -XO ) / 

(TABLE-INPUT ( I , 2 ) -YO ) ) 

IF  LIST-ALFA(I)  >0.0  AND  XO  <  TABLE-INPUT (1,1)  THEN 
LIST-AZ (I )4-LIST-ALFA(I) 

ELSE  IF  LIST-ALFA(I)  <0.0  AND  XO  >  TABLE-INPUT  (1 . 1 )  THEN 
LIST-AZ  (I)*-  LIST-ALFA  ( I )  +2 . 0*PI 
ELSE 

LIST-AZ(I)  *-LIST-ALFA(I)+PI 
END  IF 
END  IF 
END  DO 

END  INITIAL-AZIMUTHS (LIST-AZ) 

MODULE  38 

ALGORITHM  MATRIX-L (TABLE- INPUT (1,3), LIST-AZ , N ) 

ALGORITHM  ZERO(LIST-L) 

DO  FOR  I*-l  TO  10 
LIST-L(I)  4-0.0 
END  DO 
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END  ZERO(LIST-L) 

DO  FOR  J  *•  1  TO  N 
J  4-1  +  1 

LIST-L ( I ) +— TABLE-INPUT (1,3) +LIST-AZ (I )-LIST-AZ ( J ) 

END  DO 

END  MATRIX-L (LIST-L) 

MODULE  39 

ALGORITHM  SQUARED-DISTANCES (TABLE-INPUT ,XO ,YO) 

DO  FOR  I**I  TO  M 

LIST-SO ( I )  ^-(TABLE-INPUT (1,1) -XO )**2  +  (TABLE-INPUT (1,2) 
YO)**2 

END  DO 

END  SQUARED-DISTANCES (LIST-SO) 

MODULE  40 

ALGORITHM  MATRIX-A (N , TABLE-INPUT , XO , YO , LIST-SO ) 

ALGORITHM  ZERO(TABLE-A) 

DO  FOR  14-1  TO  10 
DO  FOR  J*-l  TO  2 
TABLE-A(I,J)4-0.0 
END  DO 
END  DO 

END  ZERO(TABLE-A) 

DO  FOR  I*-l  TO  N 
J+-I+1 

TABLE-A(1 ,1)  *-(YO-TABLE-INPUT ( J , 2) ) /LIST-SO (J )- 
( YO-TABLE-INPUT (1,2)) /LIST-SO ( I ) 

END  DO 

DO  FOR  I*-l  TO  N 
J  4r*  1  +  1 

TABLE-A( 1 ,2)4-  (YO-TABLE-INPUT (1,1)) /LIST-SO ( I ) - 
(XO-TABLE-INPUT (J , 1 ) ) /LIST-SO ( J ) 

END  DO 

END  MATRIX-A (TABLE-A) 


MODULE  50 

ALGORITHM  FIX-BY-TWO-RANGES -AND-ONE-AZIMUTH 
N  *-3 

PI  *-3.141592653589793 
DO  FOR  I*-l  TO  3 

INPUT  TABLE-INPUT (1,1), TABLE-INPUT (1,2), TABLE-INPUT (1,4) 
END  DO 

DO  FOR  I*-l  TO  2 

OUTPUT’ ST#’ , I, 'EAST-' , TABLE-INPUT (1,1), ’ NORT® ' , TABLE- 
INPUT  ( I , 2 ) , ’ ST  ERROR® ' , TABLE-INPUT (1,4),' METERS ’ 

END  DO 

OUTPUT ' ST#3EAST= ' , TABLE-INPUT (3,1),’ NORT® ’ , TABLE- 

INPUT  (  3 , 2  )  ,  f  ST  ERROR® ’ , TABLE-INPUT (3,4),’ DEGREES ’ 
ALGORITHM  CONVERSION-DEGREES -RADIANS (TABLE-INPUT (3,4)) 
TABLE-INPUT (3,4)*- TABLE-INPUT ( 3 , 4 )* (PI / 180 . 0 ) 

END  C0NVERSI0N-DEGREES-RADIANS(TABLE-INPUT(3 ,4) ) 

INPUT  TABLE-INPUT (1,3), TABLE-INPUT (2,3), TABLE-INPUT (3,3) 

DO  WHILE  TABLE-INPUT (1,3) ^0.0 

OUTPUT 'OBSERVED  RANGE  DISTANCES  AND  AZIMUTH  ANGLE' 

OUTPUT 'Rl=' , TABLE-INPUT (1,3), 'METERS' 

OUTPUT ' R2=' , TABLE-INPUT (2, 3) . 'METERS’ 

OUTPUT ' A= ' , TABLE-INPUT ( 3 , 3 ) , f  DEGREES ' 

ALGORITHM  CONVERSION-DEGREES-RADIANS (TABLE-INPUT (3,3)) 
TABLE-INPUT (3,3)*- TABLE-INPUT ( 3 , 3 )* (PI / 180 . 0 ) 

END  CONVERSION-DEGREES-RADIANS (TABLE-INPUT (3 ,3) ) 

MODULES  51,52,5,53,7 

INPUT  TABLE-INPUT (1,3), TABLE-INPUT (2,3), TABLE-INPUT (3,3) 
END  DO 

END  FIX-BY-TWO-RANGES-AND-ONE-AZIMUTH 
MODULE  51 

ALGORITHM  FIRST-INITIAL-POINT (TABLE-INPUT) 

E  *- TABLE-INPUT ( 1 , 3 )**2-TABLE-INPUT ( 2 , 3 )**2+ 

TABLE-INPUT ( 2 , 2 )**2-TABLE-INPUT (1,2 )**2+ 

TABLE-INPUT ( 2 , 1 )**2-TABLE-INPUT( 1 , 1 )**2 
IF  TABLE-INPUT (2,1) f TABLE-INPUT (1,1)  THEN 
El  *- ( (TABLE-INPUT (1,2) -TABLE-INPUT (2,2))/ 

(TABLE-INPUT ( 2 ,1 ) -TABLE-INPUT (1,1;) )**2+l . 0 
E2  «- (E* (TABLE-INPUT ( 1 , 2 )-TABLE-INPUT (2,2)))/ 

( (TABLE-INPUT (2,1) -TABLE-INPUT ( 1 , 1 ) )**2 )- 
2 . 0*TABLE-INPUT ( 1 , l )* ( (TABLE-INPUT (1,2) -TABLE- 
INPUT  (2,2))/ (TABLE-INPUT ( 2 , 1 )-TABLE-INPUT ( 1 , 1 ) ) )-2 . 0* 
TABLE-INPUT (1 ,2) 

E3  4r*  ( E/  (  2 . 0* (TABLE-INPUT  (2,1) -TABLE-INPUT  (1,1)))  )**2- 
( E*TABLE-INPUT (1,1))  / (TABLE-INPUT (2,1) -TABLE- 
INPUT  (1,1 ) )-TABLE-INPUT (1,3 )**2+TABLE-INPUT( 1 , 1 )**2+ 
TABLE-INPUT (1,2 )**2 


E4«-E2**2-4.0*E1*E3 
IF  E4  <  0.0  THEN 

XO  «- TABLE-INPUT (1,1) +TABLE-INPUT ( 1 , 3 )* (TABLE-INPUT (2,1) 
TABLE-INPUT (1,1))/ SQRT ( (TABLE-INPUT (2,1)- 
T ABLE-INPUT ( 1 , 1 ) **  2 + ( TABLE-INPUT (2,2)- 
TABLE-INPUT( 1 , 2 ) )**2 ) 

YO  4- TABLE-INPUT (1,2) +TABLE-INPUT (1,3)* (TABLE-INPUT (2,2) 
TABLE-INPUT (1,2))/ SQRT ( (TABLE-INPUT (2,1)- 
TABLE- INPUT (1,1) )** 2+ (TABLE-INPUT (2,2)- 
TABLE-INPUT( 1 , 2 ) )**2 ) 

ELSE  E4=0 .0  THEN 
YO*-  -E2/(2.0*E1) 

XO  *- E/ ( 2 . 0* (TABLE-INPUT (2,1 ) -TABLE-INPUT (1,1)))  + 

YO* (TABLE-INPUT (1,2) -TABLE-INPUT (2,2))/ 

(TABLE-INPUT ( 2 , 1 ) -TABLE-INPUT (1,1)) 

ELSE 

YOl*-(-E2+SQRT(E4))/2.0*El) 

X01  VE/  (2 . 0*(TABLE-INPUT ( 2 , 1 ) -TABLE-INPUT (1 ,  ’.)))  + 

Y01* (TABLE-INPUT (1,2) -TABLE-INPUT (2,2))/ 
(TABLE-INPUT(2 , 1 ) -TABLE-INPUT ( 1 , 1 ) ) 

Y02  *- ( -E2-SQRT ( E4 ) ) / ( 2 . 0*E1 ) 

X02  «-E/ (2 .0* (TABLE-INPUT (2 , 1 ) -TABLE-INPUT (1 , 1) ) )  + 

Y02* (TABLE-INPUT (1,2) -TABLE-INPUT (2,2))/ 
(TABLE-INPUT (2,1) -TABLE-INPUT (1,1)) 

IF  TABLE-INPUT ( 3, 1)»X01  AND  TABLE-INPUT (3,2) =Y01  THEN 
XO  ♦-  X02 
YO  *-Y02 

ELSE  IF  TABLE-INPUT ( 3, 1)=X0 2  AND  TABLE-INPUT (3,2)= 

Y02  THEN 
XO+-X01 
YO  *-YOl 
ELSE 

CALL  CRITERIUM (TABLE-INPUT (3,1), TABLE-INPUT (3,2), 

XOl ,Y01 , A301 ) 

CALL  CRITERIUM ( TABLE-INPUT (3,1), TABLE-INPUT (3,2), 

X02  Y02  A3Q2) 

IF  A3Q1=TABLE-INPUT(3 ,3)  AND  A30l=A302  THEN 

OUTPUT 'SOLUTION  UNDETERMINED  FOR  THAT  DATA  SET' 
PICK  UP  ANOTHER  DATA  SET 
ELSE  IF  ( (A30l-TABLE-INPUT(3 , 3 ) )**2)= 

((A302-TABLE-INPUT(3,3))**2)  THEN 
XO  *-  ( XO  1 +X0  2 )  /  2 . 0 
YO*"(Y01+Y02)/2.0 

ELSE  IF  ( (A30 1 -TABLE-INPUT (3 , 3 ) )**2 )  > 

( (A30 2-TABLE- INPUT ( 3 , 3 ) )**2 )  THEN 
XO  *-X02 
YO  *•  Y02 
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ELSE  IF  (<A30 1-TABLE-INPUT (3, 3 ))**2)  < 

( ( A302-TABLE-INPUT ( 3 , 3 ) )**2 )  THEN 
XO— XOl 
YO  — YOl 
END  IF 
END  IF 
END  IF 
ELSE 

YO  «- E / ( 2 . 0* (TABLE-INPUT (2,2) -TABLE-INPUT (1,2))) 

F  4-  TABLE-INPUT ( 1 , 3 ) ** 2 - ( TABLE-INPUT ( 1 , 2 ) -YO )**2 
IF  F  $  0.0  THEN 

XO  —TABLE-INPUT  (1,1) 

ELSE 

X01  ^TABLE-INPUT  (1 , 1)+SQRT(F) 

Y01  — YO 

X02 4-TABLE-INPUT (1,1) -SQRT ( F ) 

Y02  — YO 

IF  TABLE-INPUT (3,1) SX0 1  AND  TABLE-INPUT (3, 2 )=Y01  THEN 
X0-X02 

ELSE  IF  TABLE-INPUT ( 3, 1)=X02  AND  TABLE-INPUT (3,2)= 

Y02  THEN 
XO—XOl 
ELSE 

CALL  CRITERIUM(TABLE-INPUT (3,1), TABLE-INPUT (3,2), 

XOl , Y01 ,A30 1 ) 

CALL  CRITERIUM(TABLE-INPUT (3,1), TABLE-INPUT (3,2), 

X02  Y02  A302) 

IF  A3oi-TABLE-INPUT (3 , 3 )  AND  A30l=A302  THEN 

OUTPUT • SOLUTION  IS  UNDETERMINED  FOR  THAT  DATA  SET 
PICK  UP  ANOTHER  DATA  SET 
ELSE  IF  ( ( A3  0 1 -TABLE- IN  PUT (3,3))**2)  = 

( (A302-TABLE-INPUT ( 3 , 3 ) )**2 )  THEN 
X0*-X01 

ELSE  IF  ((A30 1-TABLE-INPUT (3, 3 ))**2)  > 

( (A302-TABLE-INPUT(3 ,3) )**2)THEN 
XO  XO  2 

ELSE  IF  ((A30 1-TABLE-INPUT (3, 3 ))**2)  < 

( (A302-TABLE-INPUT(3 , 3 ) )**2 )  THEN 
XO-  XOl 
END  IF 
END  IF 
END  IF 
END  IF 

END  FIRST  INITIAL-POINT (XO.YO) 

MODULE  52 

ALGORITHM  ITERATIONS  (TOLERANCE) 

DO  UNTIL  TOLERANCE  < 1 . 0 

MODULES  54,55,56,57,58,2,59,12,13,14 
END  DO 

END  ITERATIONS (XO,YO) 
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MODULE  53 

ALGORITHM  PRECISION (TABLE-A , TABLE-WEIGHT , TABLE-Q , LIST-L , DELTA-X , 

DELTA  Y,N,S30) 

MODULES  18,19,60,21,22 
END  PRECISION  (SU , SX, SY , SXY ,RO) 

MODULE  54 

ALGORITHM  A30( TABLE-INPUT, XO,YO) 

CALL  CRITERIUM (TABLE-INPUT (3,1), TABLE-INPUT ( 3 , 2 ) , XO , YO , A30 ) 
END  A30(A30) 

MODULE  55 

ALGORITHM  DISTANCES (TABLE-INPUT ,XO ,YO) 

S 10  «“SQRT ( (TABLE-INPUT ( 1 , 1 ) -XO ) **2+ (TABLE-INPUT (1,2)- 
Y0)**2) 

S20  *-SQRT ( (TABLE-INPUT (2,1) -XO )**2+ (TABLE-INPUT (2,2)- 
Y0)**2) 

S30*-SQRT( (TABLE-INPUT(3 , l)-XO)**2+(TABLE-INPUT(3,2)- 
Y0)**2) 

END  DISTANCES 

MODULE  56 

ALGORITHM  MATRIX-A (TABLE-INPUT , XO , YO , S 10 , S20 , S30 , A30 ) 

TABLE-A ( 1 , 1 ) — ( XO-TABLE-INPUT ( 1 , 1 ) ) / S 1 0 
TABLE-A ( 1 , 2 ) — ( YO-TABLE-INPUT ( 1 , 2 ) ) / S 1 0 
TABLE-A ( 2 , 1 )  *-  ( XO-TABLE-INPUT ( 2 , 1 ) ) / S  20 
TABLE-A (2,2)  *-( YO-TABLE-INPUT ( 2 , 2 ) ) / S 20 
TABLE-A (3,1 ) a- COSINE (A3Q-TABLE-INPUT ( 3 , 3 ) ) * ( YO- 

TABLE-INPUT (3,2)) /S30+SINE(A30-TABLE-INPUT(3 ,3) )* 
(X0-TABLE-INPUT(3, 1) )/S30 

TABLE-A (3 , 2) *-COSINE (A30-TABLE-INPUT (3,3) )* (TABLE- 
INPUT  (3  ,  1 )-X0) /S30+SINE(A30-TABLE-INPUT(3 ,3) )* (YO- 
TABLE-INPUT  ( 3 , 2 ) ) ) /S30 
END  MATRIX-A (TABLE-A) 

MODULE  57 

ALGORITHM  LIST-L (TABLE- INPUT , S10 , S20 , S30 , A30) 

LIST-L ( 1 )<-  TABLE-INPUT ( 1 , 3 ) -S 10 
LIST-L (  2  )  —TABLE-INPUT  (  2 , 3 )  -S 20 
LIST-L(3)  *-SINE(TABLE-INPUT(3 ,3)“A30)*S30 
END  LIST-L (LIoT-L) 

MODULE  58 

ALGORITHM  BEFORE-WEIGHT-MATRIX (TABLE-INPUT ( 3 , 4 ) , S30 ) 

TABLE-INPUT (3,4)**  S30*SINE(TABLE-INPUT(3 , 4) ) 

END  BEFORE-WEIGHT-MATRIX(TABLE-INPUT(3 ,4) ) 
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MODULE  59 

ALGORITHM  AFTER-WEIGHT-MATRIX (TABLE-INPUT (3, 4) ,S30) 
TABLE-INPUT (3,4) ARC  S INE (TABLE-INPUT ( 3 , 4 ) / S30 ) 

END  AFTER-WEIGHT-MATRIX (TABLE-INPUT (3 ,4) ) 

MODULE  60 

ALGORITHM  ST-DEVIATION-OF-SACH-OBS (SU , TABLE-WEIGHT , S30) 
OUTPUT 'PRECISION  OF  OBSERVATIONS' 

DO  FOR  I«-l  TO  2 

S  — ( SU / SQRT (TABLE-WEIGHT (I , J ) ) ) 

OUTPUT 'ST  DEVIATION  OF  OBS ', I S , 'METERS ' 

END  DO 

S  -ARC  SINE(SU/(SQRT(TABLE-WEIGHT(3,3))*S30))*(180.0/PI) 
OUTPUT 'ST  DEVIATION  OF  OBS  3= S DEGREES ' 

END  ST-DEVIATION-OF-EACH-OBS 

MODULE  70 

SUBROUTINE  CRITERIUM ( XS , YS , XP , YP , ASP ) 

PI-3. 14159  26535  89793 
IF  YP-YS  AND  XP  >  XS  THEN 
ASP*- PI/2.0 

ELSE  IF  YP=YS  AND  XP  <  XS  THEN 
ASP-3. 0*PI/2.0 
ELSE 

ALFA -ARC  TANGENT  ((XP-XS)  /YP-YS)) 

IF  ALFA  *0.0  AND  XP»XS  THEN 
ASP  —  ALFA 

ELSE  IF  ALFA  <0.0  AND  XP  <  XS  THEN 
ASP  —  ALFA+2 . 0*PI 
ELSE 

ASP— ALFA+PI 
END  IF 
END  IF 
RETURN 

END  CRITERIUM 
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NUMBER  OF  STATIONS*  3 


ST# 

1 

EAST* 

595794.50 

NORT* 

4055042.70 

ST 

ERROR* 

0.0200 

ST# 

2 

EAST* 

597967.80 

NORT* 

4053453.20 

ST 

ERROR* 

0.0240 

ST# 

3 

EAST* 

603425.20 

NOR  T* 

4053917.20 

ST 

ERROR* 

0.0180 

OBSERVED  AZIMUTHS 

AZIMUTH  FRCM  STATION#  1  =  76.017  DEGREES 
AZIMUTH  FRCM  STATION#  2  *  45.541  DEGREES 
AZIMUTH  FRCM  STATION#  3  =313.005  DEGREES 


ADJUSTED  COORDINATES  X*  600868.306  Y=  4056302.781 


PRECISION  OF  OBSERVATIONS 
ST  CEVIATICN  OF  OBS  1  *0.031  DEGREES 

ST  DEVIATION  OF  CBS  2  *0.038  OEGREES 

ST  DEVIATION  OF  03S  3  *0.023  OEGREES 

SX*  2.13  SY*  1.70  S  XY=  -0.862 
CORRELATION  COEFFICIENT  RG=-.24 


ERROR  ELIPSE  SEMI-AXIS  AND  ORIENTATION 
SEMI-MAJOR  AXIS  SA=  2.283 
SEMI-MINOR  AXIS  SB*  1.636 

ANGLE  FRCM  X-AXIS  TO  SA  ANT ICLOCKW I S E*157 . ODEG 


NUMBER 

CF  SEXTANT  ANGLES®  3 

ST* 

1 

EAST* 

603425.20  NORT= 

4053917.20 

ST 

ERROR= 

ST* 

2 

EAST  = 

600372.00  NORT= 

4051216.90 

ST 

ERROR* 

ST* 

3 

EAST  * 

597967.80  NORT= 

4053453.20 

ST 

ERROR* 

ST# 

4 

EAST  = 

595794.50  NOR  T= 

4055042.70 

OBSERVED 

SEXTANT 

SEXTANT  ANGLES 

ANGLE  BETWEEN 

ST# 

1 

AND 

ST# 

2  = 

49.927 

DEGREES 

SEXTANT 

ANGLE 

BETWEEN 

ST# 

2 

AND 

ST# 

3  = 

38.130 

DEGREES 

SEXTANT 

ANGLE 

BETWEEN 

ST# 

3 

AND 

ST# 

4  = 

30.396 

DEGREES 

ADJUSTED  COORDINATES  X=  600864.586  Y=  4056512.323 


PRECISION  OF  OBSERVATIONS 
ST  DEVIATION  OF  OBS  1  =0.007  DEGREES 

ST  DEVIATION  CF  GBS  2  =0.007  DEGREES 

ST  DEVIATION  CF  OBS  3  =0.007  DEGREES 

SX=  1.02  SY=  0.48  SXY=  -0.097 

CORRELATION  COEFFICIENT  R0=-.20 


ERROR  ELIPSE  SEMI-AXIS  AND  ORIENTATION 
SEMI-MAJCP  AXIS  SA=  1.050 
SEMI-MINOR  AXIS  S3=  0.480 

ANGLE  FRCM  X-AXIS  TO  SA  ANT ICLOCKW I S E=173 . 30EG 


1.0000 

1.0000 

1.0000 
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ST# 

1 

EAST* 

595794.50 

NORT* 

4055042.70 

ST 

ERROR* 

10.000 

ST# 

2 

EAST* 

603425.20 

NORT* 

4053917.20 

ST 

ERROR* 

10.000 

ST# 

3 

EAST* 

597967.80 

NOR  T* 

4053453.20 

ST 

ERROR* 

0.024 

OBSERVED  RANGE 
Rl*  5233.00 
R2*  3515. CO 

A=  45.54 


DISTANCES  AND 
METERS 
METERS 
DEGREES 


AZIMUTH  ANGLE 


ADJUSTED  COORDINATES  X*  600872.166  Y  =  4056304.120 


PRECISION  CF  OBSERVATIONS 
ST  DEVIATION  CF  08S  1  =  3.45  METERS 

ST  DEVIATION  CF  OBS  2  =  3.45  METERS 

ST  DEVIATION  OF  CBS  3  =  0.008  DEGREES 

SX=  2.86  SY=  2.88  SXY*  7.911 
CORRELATION  COEFFICIENT  RC  =  0.96 


ERROR  ELIFSE  SEMI-AXIS  AND  ORIENTATION 
SEMI-MAJOP  AXIS  SA*14.265 
SEMI-MINOR  AXIS  SB*  2.051 

ANGLE  FROM  X-AXIS  TO  SA  ANTICLOCKWISE*  45  .2DEG 


METERS 

METERS 

DEGREES 
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